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OPERATION  ANALYSIS 


INTRODUCTION 


This  report  documents  work  don*  on  the  Electro-Mechanical-Constant 
Recoil  study  performed  under  DH96  and  AH78  at  Rock  Island  Arsenal,  Rock 
Ialand,  IL,  61201. 

CONCEPT 

A self-powered,  nigh  impulse  generator  was  to  be  developed  to  interface 
with  an  analog  logic  system  to  enable  test  data  be  obtained  with  respect  to 
a constant  recoil  system.  The  telescoped  cased  consolidated  30MM  AMCAWS 
round  was  selected  due  to  its  availability  and  high  impulse  levels  (150 
lb/sec.)  Internal  reports  relating  to  constant  recoil  (firing  out  of  bat- 
tery) were  the  starting  point  of  this  program.  A self-powered  version  was 
decided  upon  since  an  externally  powered  weapon  1b  being  used  for  a recoil 
force  reduction  program  underway  at  Minneapolis-Honeywell  of  St.  Louis  Park, 
Minnesota. 


PROGRESS  SUMMARY 

The  major  objective  of  this  effort  was  to  design  a self-powered  firing 
fixture  capable  of  handling  medium  to  high  Impulses  at  varying  rates  of 
fire.  The  initial  criteria  was  425  SPM  and  150  lb/sec  Impulse.  The  "ideal" 
minimum  recoil  force  is  obtained  by  accelerating  the  recoiling  masses  in  the 
counter  recoil  (or  forward)  direction,  up  to  a velocity  that  will  cancel 
one-half  of  the  round  Impulse  of  150  lb/sec.  The  recoil/counterrecoil  stroke  of 
the  main  mass  is  5 Inches  while  the  slide  travels  an  additional  9 1/2  inches 
or  a total  of  14  1/2  Inches.  This  extra  slide  travel  is  for  feeding  clear- 
ances and  chamber  locking.  An  analysis  of  the  operating  principle  is  in 
the  Appendix. 


TEST  DATA 


Single  round  firings  utilizing  a set  of  AMCAWS  30  recoil  buffers  were 
performed  to  determine  proper  gas  orifice  size  in  the  barrel  which  is  .156 
in.  and  required  slide  spring  energies  to  obtain  reliable  sear-up,  which  turned 
out  to  be  3500  in/lb.  The  slide  recoil  stroke  of  9 1/2  in.  took  33ms 
which  is  good.  The  counterrecoil  stroke  of  the  slide  required  95ms  which 
is  too  slow.  Increasing  spring  energy  further  is  detrimental  to  the  recoil 
stroke  (See  Problem  Solutions  #7). 

Results  of  18  Single  Round  Firings. 

1.  Firing  pin  solenoid  is  inadequate  due  to  low  energy  and  high  pull-in 
time  of  73ms  plus  drop-out  time  of  137ms.  The  cycle  time  of  the  solenoid 
alone  is  210ms  = 286  cyc/min. 

2.  Bolt  extension  and  stud  rollers  were  too  soft. 

3.  Feed  space  between  bolt  and  barrel  face  was  too  small. 

4.  Firing  pin  assembly  was  too  heavy  Increasing  action  time. 

5.  Piston  configuration  is  satisfactory. 

6.  Gas  port  of  .156  dia.  (3  pics)  in  barrel  at  44"  is  satisfactory. 

7.  Drive  springs  for  slide  require  3500  in/lb  energy. 

8.  Baseplate  must  be  stiffened  to  assist  sear  operation. 

9.  Gun  cycle  time  is: 


Recoil  stroke 

35ms 

9 

Feed  time 

30ms 

Counterrecoil  stroke 

95ms 

Cycle  time 

160ms 

A firing  rate  of  375  SPM  results. 
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10.  The  firing  solenoid  + mechanism  + round  action  times  must  he  concur- 
rent with  the  counterrecoil  stroke  time  which  creates  logic  lead-time  prob- 
lems . 

11.  The  main  sear  required  excessive  force  to  release  it  requiring  additional 
analysis . 

« 12.  Installation  and  removal  of  the  slide  springs  is  extremely  difficult. 

13.  Main  mass  charger  pinion  skips  teeth  due  to  rack  bowing  under  the  1200 
lb.  main  mass  buffers. 

14.  Loss  of  head  space  after  each  firing  occurred  due  to  lack  of  adequate 
barrel  to  receiver  locking  method. 

15.  Alignment  of  gas  cylinder,  to  center  guide,  to  piston,  to  slide,  is 
critical. 

PROBLEM  SOLUTIONS 
(Refer  to  test  results  above) 

1.  Solenoid  and  firing  mechanism  reworking  reduced  pull-in  time  to  20ms  and 
drop-out  time  to  100ms  for  cycle  time  of  120ms  which  was  still  too  long, 
therefore  a redesigned  firing  pin  mechanism  and  actuating  device  were  designed 
and  fabricated.  See  drawing  76F40268  for  new  design.  Parts  have  been  fab- 
ricated ana  assembled. 

2.  The  bolt  extension  and  stud  rollers  were  hardened. 

3.  Space  was  increased  by  adding  \ inch  to  recoil  motion  of  slide  there 
by  having  a 9 3/4  inch  slide  motion  relative  to  the  receiver. 

m 

4.  See  1 above. 

5.  N/A 

6.  N/A 

7.  Drive  springs  for  the  slide  are  difficult  to  Install  and  remove  due  to 
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the  large  forces  Involved.  Since  the  reduction  of  spring  load  increases 
slide  travel  time  which  is  presently  too  great  an  assembly  procedure  to 
offset  this  problem  was  devised.  See  II  E of  assembly  procedures. 

8.  A new  Firing  Table  was  fabricated  that  incorporated  all  the  functional 
needs  of  the  original  base  plate. 

9.  N/A 

10.  See  1 above. 

11.  A new  main -sear  mourted  in  the  right  rear  support  plate  is  actuated 
by  an  air  cylinder. 

12.  See  7.  above. 

13.  The  pinion  and  rack  no  longer  charge  against  the  1200  lb  main  mass 
buffer.  A separate  main  charger  was  designed. 

14.  A positive  locking  device  must  still  be  designed. 

15.  No  analysis  has  been  undertaken  on  this  problem. 
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AMC  30  ACT  ASSEMBLY  PROCEDURE 


I.  Firing  Table  Components  Assembly 

A.  Crank  Sub-Assembly  (73B40299) 

1.  Pin  crank,  handle,  bolt,  and  pinion  shaft  (73B40235)  together. 

2.  Locate  latch  (75B40188)  (large  opening  up)  in  lock  (75B40189) 
and  insert  pinion  shaft. 

3.  Place  pinion  (73C40242)  over  shaft  and  pin  in  place. 

B.  Slide  Charger  Sub-Assy 

1.  Locate  charger  block  (75C40187)  on  table  and  secure  with  four 
bolts. 

2.  Insert  spring  (73B40238)  in  block. 

3.  Compress  spring  with  shaft  of  crank  sub-assy  and  fasten  in  place 
with  two  bolts.  (See  A for  Crank  Sub-Assy) 

C.  Rear  Support  Assembly 

1.  Locate  brace  (75C40182)  on  firing  table  and  bolt  in  position. 

2.  Assemble  left  and  right  rear  supports  to  brace.  Assure  6.00 
inch  space  between  rail  guides. 

3.  Assemble  pivot  block  (75B40119),  main  sear  (75C40120),  leaf 
spring  (75B40118),  and  retaining  plate  (75C40121),  in  that 
order  to  right  rear  support. 

4.  Holding  main  sear  retracted,  locate  stud  (76B40232)  and  insert 
cross  pin  (76B40230).  Allow  sear  to  return  into  sideplate  which 
retains  pin  in  position. 

AFTER  II  COMPLETE  THE  FOLLOWING  STEPS  5 & 6 

5.  With  spacer  plate  in  position,  screw  air  cylinder  on  stud  (4 
above)  until  the  main  sear  protrudes  into  plate  mid-groove 
1/8  inch. 

6.  Continue  rotation  until  air  inlet  is  in  up  position,  then  se- 
cure cylinder  to  side  plate  with  2 bolts. 


D.  Front  Support  Assy 

1.  Locate  front  spacer  (7GR40229),  support  (72040225),  and  straps 
(76B40231)  in  that  order  on  table  and  bolt  down. 

Weapon  Assembly 

A.  Chamber  - Bolt  - Firing  Pin  Sub  Assy  (76D40270) 

1.  Start  with  chamber. 

2.  Insert  bolt  (75B40238)  into  chamber. 

3.  Assemble  firing  pin  (76A40240),  tube  (76B40239) , and  retaining 
pin. 

4.  Insert  spring  (75A40183)  and  then  3 above  into  bolt. 

5.  Place  backplate  sub-assy  (73C40297)  on  bolt. 

6.  Insert  bolt  extension  (76C40237)  in  backplate. 

7.  Install  crosspin  (73B40183)  thru  bolt  and  tube  holding  it  in 
position  with  spring  (73B40181-C) . 

8.  Compress  spring  with  base  (73B40186)  and  insert  retaining  pin 
thru  tube  and  plug. 

B.  Barrel  Sub-Assy 

1.  Locate  headspace  nut  (73B40138)  on  barrel. 

2.  Place  central  guide  (73C40208)  and  nut  (73B40209)  on  barrel 
(DO  NOT  TIGHTEN) . 

3.  Screw  plug  (73B40212)  with  seal  (73B40207)  into  gas  cylinder. 

4.  Slip  gas  cylinder  and  nut  (73B40211)  onto  barrel  (DO  NOT  TIGHTEN). 

5.  Insert  gas  piston  (73D40147)  into  cylinder. 

C.  Barrel-Receiver  Group  Sub-Assy 

1.  Install  receiver  upside  down  in  rear  support  (Stud  slot  up). 

2.  Screw  barrel  3 inches  into  receiver. 

3.  Install  rear  slide  guide  (73C40144). 

4.  Insert  chamber  sub-assy  into  receiver  screwing  in  backplate 
until  it  is  flush  with  the  rear  of  the  receiver. 


5.  Screw  in  the  backplate  alignment  screw  (73B40177). 

6.  Rotate  bolt  so  longitudinal  slot  lines  up  with  the  receiver 
stud  slot. 

7.  Screw  on  housing  (75D40138)  hand  tight. 

8.  Install  stud  (73B40159)  and  rollers  (73B40158)  (thin  roller 
first)  . 

9.  Assemble  round  stop  pawl  (73B40161). 

10.  Insert  headspace  gage  and  lock  chamber  in  forward  position  (gage 
set  at  6.305). 

11.  Screw  in  barrel  tightly  and  secure  with  head  space  nut. 

12.  Remove  headspace  gage. 

13.  Orientate  central  guide  with  holes  up. 

14.  Slip  on  two  spring  housings  (73C40215) . 

15.  Lay  slide  sub-assy  (See  D Below)  against  receiver  locating  cam 
path  over  the  stud  roller. 

16.  Connect  slide  to  rear  slide  guide  and  secure  with  pip  pin. 

17.  Connect  slide  to  piston  with  retaining  pin. 

D.  Slide  Sub-Assembly 

1.  Assemble  slide  safety  using  lock  (75C40162),  torsion  spring 
(75B40161),  and  shoulder  stud  (75B40160)  to  the  slide. 

2.  Insert  pip  pin  thru  sear  into  slide. 

3.  Install  buffer  lug  (73B40142)  and  spring  and  pin  in  place. 

4.  Install  sear  assembly  (73B40298)  in  the  slide. 

5.  Locate  sear  holder  (77C40003)  on  slide  to  compress  sear  assembly. 

E.  Slide  Spring  Installation 

1.  Turn  barrel-receiver  assy  over  in  rear  supports. 

2.  Assemble  receiver  strap  (77B40001)  to  back  of  rear  supports. 

3.  Clamp  slide  support  (77C40002)  2"  forward  of  rear  support. 


4.  Engage  pinion  with  slide  rack  and  crank  slide  forward  until 
rear  slide  guide  pip  pin  is  accessible. 

5.  Disengage  rear  slide  guide  and  while  holding  receiver  against  2. 
above  continue  to  crank  the  slide  forward  until  the  pinion  and 
rack  are  almost  disengaged. 

6.  Remove  slide  safety  pip  pin. 

7.  Assemble  drive  springs  on  guides  (73C40192)  and  insert  into 
spring  housings. 

8.  Insert  two  pins  (77B40004),  in  charger  block. 

9.  Crank  slide  assembly  back  until  slide  safety  (77C40002)  can  be 
installed.  Install  it. 

10.  Engage  rear  slide  guide  to  slide  and  repin. 

11.  Locate  the  rail  slide  guide  (73B40213)  and  fasten  in  position. 

12.  Disengage  slide  safety  and  pin  in  place. 

13.  Remove  slide  support. 

14.  Remove  pin  in  buffer  lug. 

15.  Remove  sear  holder. 

16.  Remove  two  pins  in  charger  block. 

17.  Felieve  load  on  slide  safety  in  8.  above  and  remove  it. 

CAUTION:  Crank  will  have  a spring  load  which  must  be  relieved 

gradually. 

18.  Align  central  guide  and  gas  cylinder,  securing  them  in  place 
with  their  respective  nuts. 

Firing  Pin  Actuator  Sub  Assy  (76F40268) 

1.  Slip  solenoid  mount  (75C40137)  over  housing  (75D40138). 

2.  Align  holes  and  insert  tow  studs  (76B40241)  finger  tight. 

3.  Insert  two  pins  (75A40132)  into  guide  (75B40136). 

Install  striker  (75B40135,  spring  (75B40134)  and  guide  into 
housing  and  cam  into  captive  position. 


4. 


5.  Rotate  housing  so  two  studs  are  up  and  install  screw  (MS51962-52) 
in  tapped  hole  between  the  two  studs. 

6.  Place  shim  (76B40245),  plunger  (75C40133),  spring  (see  11701134) 
and  solenoid  into  mount  and  secure  with  nut  (MS  35649-202)  and 
cap  screw  (MS16995-36) . 

7.  Rotate  housing  180°  and  install  spring  (75A40185)  on  long  end 
of  pin  (76A40233) . 

8.  Insert  pin  with  spring  down  into  hole  in  solenoid  mount. 

9.  Using  cap  screws  (MS16995-42) , assemble  cover  (76A40234),  pin 
holder  (76A40235) , pin  (76A40269),  and  block  (76A40236)  together 
In  that  order. 

10.  Locate  9.  over  pin  in  8.  above  and  secure  to  mount  with  cap 
screw. 

11.  Tighten  studs  (76Br0241)  and  firing  mechanism  to  backplate. 

Assembly  and  Calibration  of  Feeder 

1.  Screw  on  front  bracket. 

2.  Insert  rod  (73B40168). 

3.  From  front,  insert  rod  (73B40248). 

4.  Screw  on  retainer  pin  (1  of  4)  (7340166). 

5.  Slip  on  torsion  spring,  feed  arm,  feed  release. 

6.  Screw  retainer  pin  (2  9f  4)  thru  feed  arm. 

7.  Screw  retainer  pin  (3  of  4)  thru  feed  release  making  sure  front 
rod  is  captured.  Check  rotation  of  feed  arm  to  assure  that 
front  rod  is  not  180°  off. 

8.  Screw  on  rear  bracket. 

9.  Slip  in  spring  (73B40165). 

10.  Screw  retainer  pin  (4  of  4)  thru  rear  bracket,  compressing 
spring  (73B40165). 

11.  Screw  on  cartridge  guide  (to  receiver). 

12.  Screw  on  rammer  base. 


H.  Assembly  Procedure,  Magazine 

1.  Begin  with  housing  weldment. 

2.  Insert  spring  with  the  end  loop  in  the  upper  hole. 

3.  Place  the  follower  assembly  in  the  lower  chute,  against  the  stop. 

4.  Insert  the  shaft,  short  end  first. 

5.  Put  arm  on  the  shaft  and  rotate  CCW  until  spring  is  wound  2-1/2 
turns. 

6.  Insert  the  hinge  pin  (MS9848-20)  thru  the  end  holes  in  the  arm 
and  the  follower  assembly.  Add  the  cotter  pin  (inside). 

7.  Add  the  hub  assy  and  screw  on,  with  9 self-tapping  screws. 

8.  Test  by  loading  25  dummy  rounds.  Rounds  should  feed  freely. 


OPERATION  OF  FIXTURE 


Push  in  pinion  crank  to  engage  rack.  Lock  in  position  with  latch. 
Withdraw  pip  pin  from  slide  safety  lock. 

Crank  back  slide  until  slide  sear  engages  receiver. 

Disengage  charger  pinion  from  rack,  slide  sear  holds  back  slide  along 
with  safety  lock. 

Position  main  mass  charger  on  rear  support  and  crank  back  until  main 
sear  engages  the  receiver  rail.  Manually  pull  back  feeder  and  insert 
cartridge. 

Remove  main  mass  charger  and  repin  slide  safety  lock. 

Main  sear  is  released  by  activating  the  air  cylinder. 

The  slide  sear  is  released  by  being  tripped  by  a properly  positioned 
plate  under  the  slide. 
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INPUT  PARAMETERS 


1. 

PRMT(l) 

Initial  integration  time  (sec) 

2. 

P RMT ( 2 ) 

Final  integration  time  (sec) 

3. 

P RMT ( 3 ) 

Integration  time  step  (sec) 

4. 

PRMT (4) 

Allowed  error  (dimensionless) 

5. 

PRMT (5) 

Control  number  (set  to  zero) 

6. 

QRANG(l) 

Expected  displacement  range  of  receiver 
(in) 

7. 

QRANG (2) 

Expected  displacement  range  of  slide  (in) 

8. 

QDRANG(l) 

Expected  velocity  range  of  receiver 
(in/sec) 

9. 

QDRANG(2) 

Expected  velocity  range  of  slide 
(in/sec) 

10. 

Q(l) 

Initial  displacement  of  receiver  (in) 

11. 

Q (2) 

Initial  displacement  of  slide  (in) 
(if  not  in  locked  position) 

12. 

QD(1) 

Initial  velocity  of  receiver  (in/sec) 

13. 

QD(2) 

Initial  velocity  of  slide  (in/sec) 

14. 

TRANG 

Expected  temperature  range  of  cylinder 
gas  (° R) 

15. 

DELC , 6c 

Angle  chamber  is  when  it  becomes 
locked  (rad) 

16. 

Qcfl  » qCFL 

Displacement  of  chamber  when  fully 
locked  (in) 

17. 

QC I L , j ^ 

Displacement  of  chamber  when  locking 
starts  (in) 

18. 

EPSCI,  e CI 

Distance  to  represent  sharp  change  in 
slope  (in) 

J 


19. 

BCI , bCI 

Minor  axis  of  ellipse  for  "sharp 
slope"  (in) 

20. 

BCF,  b.,c 
Lr 

Minor  axis  of  ellipse  for  "sharp 
slope"  (in) 

21. 

UFL,  UpL 

(q2~qc)  when  chamber  is  fully 
locked  (in) 

22. 

UIL,  UIL 

(q2*qc)  when  chamber  is  starting 
to  lock  (in) 

23. 

EPSQF , e^£ 

Distance  to  describe  a sharp  change 
in  slope  (in) 

24. 

BQF.  bqf 

Minor  axis  of  ellipse  for  "sharp 
slope"  (in) 

25. 

ERROR,  ec 

Error  allowed  for  computing  constraint 
function  (rad) 

26. 

DELA,  6a 

Distance  slide  travels  before  gas  is 
cut  off  (in) 

27. 

DELB,  6b 

Distance  slide  travels  when  gas  flows 
back  to  barrel  (in) 

28. 

AO,  Ao 

Throat  area  (in2) 

29. 

co,  co 

Orifice  coefficient  (dimensionless) 

30. 

CVB’  Cvb 

Specific  heat  at  constant  volume  of  barrel 
gas  in 2 

Sec  2 - °R 

31. 

CPB,  Cpb 

Specific  heat  at  constant  pressure  of 
barrel  gas 

Sec2  °R 

32. 

VO,  vo 

Initial  cylinder  volume  (in3) 

33. 

ACYL,  Ac 

Cylinder  cross  sectional  area  (in2) 
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34. 

PA’  Pa 

Atmospheric  pressure  ( PS I ) 

35. 

TA,  Ta 

Atmospheric  temperature  (°R) 

36 . 

TP0RT’  Tport 

Time  for  projectile  to  reach  port  (sec) 

37. 

IBM,  Tbm 

Temperature  of  chamber  gas  when 
projectile  leaves  barrel  (°R) 

38. 

PBM*  Pbm 

Pressure  of  chamber  gas  when  projectile 
leaves  barrel  (PSI) 

39. 

NP 

# of  discrete  (time  vs  chamber  pressure) 
points 

40. 

TIMEP  (I ) 

Time  coordinates  (sec)  after  projectile 
passes  port 

41. 

PB  ( I ) 

Pressure  coordinates  (PSI)  after 
projectile  passes  port 

42. 

NT 

# of  (time  vs  chamber  temperature) 
coordinates 

43. 

TIMET (I)  . 

Time  coordinates  (sec)  after  projectile 
passes  port 

44. 

TB  (1 ) 

Temperature  coordinates  (°R)  after 
projectile  passes  port 

45. 

NCHAM 

# of  coordinates  describing  chamber 
pressure  force 

46. 

CHAMT ( I ) 

Time  coordinates  (sec) 

47. 

CHAMF(I) 

Force  coordinates  (lb) 

48. 

QiREF,  qlref 

q.  - coordinates  describing  recoil 
1force  (in) 

49. 

F1rEF,  Flref 

Force  corresponding  to  9iref  (lb) 

50. 

Q2RFF , q2ref 

- coordinates  describing  recoil 
Aforce  (in) 

51. 

F2REF,  F2ref 

Force  corresponding  to  q?re£  (lb) 
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52. 

Q3RFF , q3ref 

53. 

F3REF , F 3ref 

54. 

CR“ * cr2 

55. 

PLnF’  Fpld(f) 

56. 

Q0F’  qo(f) 

57. 

PELSF, 

% loss  (f) 

58. 

XK1F , K1(f) 

59. 

EPSF,  E(f) 

60. 

PLDR,  Fp id (r ) 

61. 

Q0R’  qo(r) 

62. 

PELSR, 

% loss  (r) 

63. 

XK1R,  K1(r) 

64. 

EPSR,  E(r) 

65. 

IDE LAY,  Tdelay 

66 . 

FSRFF*  Fsref 

67. 

XKS1 , Ksl 

q^  - coordinates  describing  recoil 
force  (in) 

Force  corresponding  to  q3ref-  (lb) 

Viscous  frictional  coeff icient (receiver) 
.lb-sec-, 

^ in  ’ 

Preload  force  on  forward  buffer  (receiver) 
(lb) 

q^  - coordinates  when  receiver  contacts 
forward  buffer  (in) 

Percent  energy  loss  from  impact  of 
forward  buffer  (%) 

Spring  rate  of  forward  buffer  (i^) 

Distance  required  for  force  to  build  up 
to  the  preload  (in) 

Preload  force  on  rear  buffer  (receiver) 
(lb) 

(-q^)  - coordinates  when  receiver 
contacts  rear  buffer  (in) 

Percent  energy  loss  from  impact  of  rear 
buffer  (%) 

Spring  rate  of  receiver  rear  buffer  (i^) 

Distance  required  for  force  to  build  up 
to  preload  of  rear  buffer  (in) 

Delay  between  firing  and  chamber  force 
starting  (sec) 

Slide  drive  spring  preload  (lb) 

Spring  rate  of  slide  drive  spring  (ij^) 
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68  . 

CS2,  Cs2 

Viscous  frictional  coefficient 
(slide)  (M) 

69. 

PLDSR’  Fpld(sr) 

Preload  force  on  slide  rear  buffer  (lb) 

70. 

Q°SR,  Qo(sr) 

( -q 2 ) - coordinate  when  slide  contacts 

rear  buffer  (in) 

71. 

PELSSR, 

% loss  (sr) 

Percent  energy  loss  from  impact  of  rear 
buffer  spring 

72. 

XK1SR , K1(sr) 

Spring  rate  of  slide  rear  buffer  spring 

(±i) 

v in' 

73. 

EPSSR,  E(sr) 

Distance  required  for  force  to  build  up 
to  preload  (i^-) 

74. 

PLDSS*  Fpld(ss) 

Preload  force  on  sear  buffer  spring  (lb) 

75. 

Q°ss,  Q0(ss) 

q2  - coordinate  when  slide  contacts 
sear  buffer  (in) 

76. 

PELSSS, 

* loSS(ss) 

Percent  energy  loss  from  impact  of 
sear  buffer 

77. 

XK1SS,  K1(ss) 

Spring  rate  of  sear  buffer 

78. 

EPSSS , E(ss) 

Distance  required  for  force  to  build  up 
to  preload  (in) 

79. 

PLDSF,  Fpld(sf) 

Preload  on  slide  forward  buffer  spring 
(lb) 

80. 

Q°SF’  Qo(sf) 

q2  - coordinate  when  slide  contacts 

forward  buffer  (in) 

81. 

PELSSF, 

% loss  (sf) 

Percent  energy  loss  from  impact  of 
forward  buffer 

82. 

XK1SF , K1(sf) 

Spring  rate  of  forward  buffer  (|^) 
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8 3. 

EPSSF*  H ( s f ) 

Distance  required  for  force  to  build  up 
to  preload  (in) 

84. 

PLDCR , *'pl<ncrj 

Preload  on  chamber  rear  buffer 

85. 

QOCR,  Qo(cr) 

(-q  ) - coordinate  when  chamber  contact: 

rear  buffer  (in) 

86. 

PELSCR, 

1 l0SS  (cr) 

Percent  energy  loss  from  impact  of  rear 
buffer  spring 

87. 

XK1CR,  K1(cr) 

Spring  rate  of  chamber  rear  buffer 
spring  (i£) 

88. 

EPSCR,  E(cr) 

Distance  required  for  force  to  build 
up  to  the  preload  (in) 

89. 

cc,  cc 

Viscous  frictional  coefficient  for 
chamber  (-1— j^?-c) 

90. 

G 

Gravitational  constant  (— — ) 

sec2 

91. 

THETA 

Initial  elevation  angle  (not  applicable 
in  model  with  vehicle  motion) (Rad) 

92. 

Q1REL 

Position  of  receiver  when  slide  is 
released 

93. 

T1 

Time  delay  (start  of  chamber  force  to 
peak  value) (sec) 

94. 

NINT 

Printing  frequency  of  forces,  pressures 
etc . 

95. 

I DELAY 

Control  # 

1 - if  correcting  for  time  delay  in 
firing  logic 

0 - if  approximating  recoil  force 
with  non-correctant  function 

96. 

IPCNT 

Control  # (1  if  approximating  recoil 
force  with  a constant  function) 
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97. 

I EC  NT 

Control  II  (1  if  correcting  for 
c 1 eva  t i on ) 

98. 

1 1CNT 

Control  H U if  correcting  for 
inertial  forces) 

99. 

I WRITE 

Control  H (1  if  desire  output  on  disk) 

100. 

I LOCK 

Control  # (1  if  chamber  starts  off 
being  locked) 

101. 

WSL 

Weight  of  slide  (lb) 

102. 

wen 

Weight  of  chamber  (lb) 

103. 

WREC 

Weight  of  receiver  (lb) 

104  . 

CIXX 

Moment  of  inertia  of  chamber  (about 
translational  axis ) ( lb- in  - sec2 ) 

105. 

Q2PTI  (I) 

q 2 ~ coordinate  when  "zth"  constant 
dissipative  force  begins  (in) 

106. 

Q2PTI1 (I) 

- coordinate  when  "zth"  constant 
dissipative  force  ends  (in) 

107. 

EPSQ21  (I) 

Build  up  distance  for  "zth"  force 
(slide) (in) 

108. 

ERICQ2 ( I ) 

Value  of  force  (slide)  (lb) 

109. 

ICNQ21  (I) 

Control  # (if  1,  equal  and  opposite 
force  acts  on  receiver) 

110. 

NQ2 

Number  of  forces  opposing  when  velocity 
is  negative 

Ill  . 

MQ2 

Number  of  forces  opposing  when  velocity 
is  positive 

112. 

Q1PTI  (I) 

q j - coordinate  when  "zth"  constant 
dissipative  force  begins  (in) 

113. 

QlPTIl(I) 

q^  - coordinate  when  "zth"  constant 
dissipative  force  ends  (in) 
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114. 

F.PSQ11  (I) 

Build  up  distance  for  "zth"  force 
(receiver) ( in) 

115. 

FRICQl(I) 

Value  of  force  (rece i ver ) (lb) 

116. 

ICNQ1I  (I) 

Control  # (no  equal  and  opposite  force 
acting)  Arbitrary 

117. 

NQ1 

Number  of  forces  opposing  motion  when 
velocity  is  negative 

118  . 

MQ1 

Number  of  forces  opposing  motion  when 
velocity  is  positive 

119. 

QCPTI (I) 

qc  - coordinate  when  "zth"  constant 
dissipative  force  begins  (in) 

120. 

QCPTI1  (I) 

qc  - coordinate  when  "zth"  constant 
dissipative  force  ends  (in) 

121. 

FPSQCI  (I ) 

Build  up  distance  for  "zth"  force  (in) 

122. 

FRICQC(I) 

Value  of  force  (chamber)  (lb) 

123. 

ICNQCI  (I) 

Control  # (1  if  equal  and  opposite 
force  acts  on  receiver) 

124. 

NQC 

Number  of  forces  opposing  motion  when 
velocity  is  negative 

125. 

MQC 

Number  of  forces  opposing  motion  when 
velocity  is  positive 

126. 

TCPTI(I) 

ec  - coordinate  when  "zth"  constant 
dissipative  force  begins  (rad) 

127. 

TCPTI1  (I) 

oc  - coordinate  when  "zth"  constant 
dissipative  force  ends  (rad) 

128. 

EPSTCI (I) 

Build  up  distance  for  "zth"  force  (rad) 
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129. 

FRICTC(I) 

Value  of  torque  (in-lb) 

130. 

1 CNTC 1(1) 

Control  # (no  equal  and  opposite  force) 
Arbitrary 

131 . 

NTC 

Number  of  forces  opposing  negative 
angular  velocity 

132. 

MTC 

Number  of  forces  opposing  positive 
angular  velocity 

133. 

UFR 

Distance  slide  is  forward  before  weapon 
can  fire 

134  . 

I STATE 

Printing  frequence  of  state  variables 

135. 

I FILE 

File  # for  data  to  be  plotted 

136. 

I FI  LED 

File  # for  CRT  if  running  from  that 
terminal  (use  6 for  printer  and  10 
for  CRT) 

137. 

IMEAS 

Set  equal  to  1 if  the  C.M.  is  used  in 
the  control  equations  (otherwise 
receiver  state  is  the  measured 
variable) 

138. 

ISPECI 

If  using  center  of  mass  as  the  measured 
variable  for  control  equation  and  we 
desire  the  state  of  the  center  of  mass 
to  not  be  returned  to  the  initial 
displacement  with  zero  velocity,  then 
set  to  one 

139. 

QCMS 

If  ISPECI  = 1,  the  center  of  mass  will 
return  to  this  displacement  after 
firing 

140. 

QDCMS 

If  ISPECI  = 1,  the  center  of  mass  will 
return  to  this  velocity  after  firing 

141  . 

NO  TAR 

Set  equal  to  one  if  no  specified  target 

142. 

XTAR 

x - coordinate  of  the  target  (m) 
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143. 

YTAR 

y - coordinate  of  the  target  (m) 

144  . 

ZTAR 

z - coordinate  of  the  target  (m) 

NOTE:  See  "FUE  Air  Defense  Evaluation" 

for  flight  path  parameters 

145. 

THETAO 

Initial  elevation  angle  (if  no  target) 
(rad) 

146. 

THETPO 

Initial  elevation  rate  (if  no  target) 
(rad/sec) 

147. 

THEPPP 

Elevation  acceleration  (if  no  target) 
(rad/sec  2 ) 

148. 

PSIO 

Initial  azimuth  angle  (if  no  target) 
(rad) 

149. 

PSIPO 

Initial  azimuth  rate  (if  no  target) 
(rad/sec) 

150. 

PSIPPP 

Azimuth  acceleration  (if  no  target) 
(rad / sec2) 

151  . 

XSP 

x - coordinate  of  slide  C.G.  relative 
to  receiver  C.G.  (initial ) (in) 

152. 

YSP 

y - coordinate  of  slide  C.G.  relative 
to  receiver  C.G.  (initial ) (in) 

153. 

ZSP 

z - coordinate  of  slide  C.G.  relative 
to  receiver  C.G.  ( in i t ial ) (in) 

154. 

XRP 

x - coordinate  of  receiver  C.G.  relative 
to  center  of  rotation  of  the  turret 
(initial)  (in) 

155. 

YRP 

y - coordinate  of  receiver  C.G.  relative 
to  center  of  rotation  of  the  turret 
(initial) (in) 

156. 

ZRP 

z - coordinate  of  receiver  C.G.  relative 

to  center  of  rotation  of  the  turret 
(initial) (in) 
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157. 

XCP 

X 

- coordinate  of  chamber  C.G.  relative 
to  receiver  C.G.  (init ial) (in) 

158. 

YCP 

y 

- coordinate  of  chamber  C.G.  relative 
to  receiver  C.G.  (initial) (in) 

159. 

ZCP 

z 

- coordinate  of  chamber  C.G.  relative 
to  receiver  C.G.  (initial) (in) 

160. 

CIYY 

Principle  moment  of  inertia  of  the 
chamber  (lb-in-sec2) 

161 . 

CIZZ 

Principle  moment  of  inertia  of  the 
chamber  (lb-in-sec2) 
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FLIGHT  PATH  DATA 


1 . NS EG 

2.  X ( 1 ) 

3.  Y(l) 

4.  Z(l) 

5.  THETA ( 1 ) 

6.  VS(1) 

7.  VZC1) 

8.  TSTART  (1) 

9.  ITYPE 

10.  AZ 

11.  S 

12.  AS 

13.  RAD 

14.  VDOT 


Number  of  segments  in  flight  path 

x coordinate  at  start  of  flight  (m) 

y coordinate  at  start  of  flight  (m) 

z coordinate  at  start  of  flight  (m) 

Heading  coordinate  at  start  of  flight 

(degrees) 

Horizontal  velocity 
Vertical  velocity 

Time  at  start  of  flight  path  (sec) 

Control  # 

Line  segment  = 1 
Circular  segment  = 2 
Spiral  segment  = 3 

Vertical  acceleration  (g) 

Segment  length  (m)  (if  ITYPE  = 1) 

Turn  angle  (degrees)  (if  ITYPE  f 1) 

Horizontal  acceleration  (g)  (if 
ITYPE  = 1) 

Radial  accelration  (g)  (if  ITYPE  f 1) 

Turn  flag  for  ITYPE  = 2 or  3 
1 * right 
-1  = left 

Tangential  accelration  (g)  (ITYPE  = 3) 
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CONTROL  EQUATION  ANALYSIS 


The  concept  of  measuring  the  displacement  and 
velocity  of  the  recoiling  parts  and  computing  a 
parameter  such  that  when  it  becomes  greater  than 
zero  the  gun  will  fire  and  return  the  recoiling 
parts  back  to  the  origin  will  be  developed  for  a 
simple  system  and  then  corrections  will  be  made 
to  the  simple  system  to  account  for  elevation 
angle,  firing  delay,  accelerating  reference  frame, 
etc. 


Typical  Force  vs.  Deflection  Curve 
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Operating  range  chosen  to  be  as  constant  as 
possible • 


Vs  - Potential  at  sear 


Vf . - Potential  when  fired 

r ire 

• 

Xf-  , Xf.  - Displacement  and  velocity  at 
tire  tire  time  of  fire 

Since  the  impulse  applied  to  the  weapon  is  a 
delta  function,  the  potential  will  not  have  changed 
at  firing  and  only  the  velocity  will  have  decreased 
by 


I_ 

M_ 


Since  the  system  is  conservative  from  the  time  after 
fire  to  sear  and  at  sear  we  desire  a zero  velocity, 
we  have  from  the  energy  equation: 


(l.a) 


'fire  * 7Mr  <Xfire  ' ST> 1 


V. 


To  determine  the  time  of  fire,  we  compute 


(2. a) 


and  fire  when 


(3. a) 


V ♦ 


i"  (i  - ir>2  ivs 
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In  the  case  of  a constant  force 


Sea 


(4. a)  V * F (x  - x)  (in  operating  range) 

From  equations  (3. a)  and  (4. a)  fire  when 


(5. a) 


-Fx  ♦ \ Mr  (x  - J-)2  < 0 


From  equation  (5. a)  we  have: 


> 1 Mr 

- 7 r 


(x  - 


sr* 

r 


For  a given  recoil  force  (F)  , impulse  (I),  and 
mass  (Mr)  the  minimum  recoil  distance  is  found  by 

finding  when  the  right  side  of  the  above  equation 
is  a minimum. 


Minimum  occurs  when 

x “ 2ST  " *fire 


and 


min 


1 

* 


Mr 

V~ 


& 2 

r 


U 


Recoil  distance  for  conventional  artillery 


fire 


0)  is : 


x 


1 

I 


M 

r 

T~ 


2 


The  minimum  is  the  "firing  out  of  battery 
principle"  and  for  a given  impulse,  mass,  and  recoil 
force,  only  requires  a 1/4  of  the  conventional  weapon 
recoil  distance. 
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If  the  spring  has  a non-zero  slope: 


From  equation  (3. a),  fire  when 
(6. a)  - | (xV  ♦ | (x  - x')2  * 7 Mr  U - i-)2  I 0 


or 

(7. a)  (x  - 2x  ) ♦ \ Mr  (x  - jjp) 2 1 0 

If  x'-*«°  and  Kx'-*-F,  then  K-»0  and  the  control 
equation  (7. a)  reduces  to: 

(8. a)  -Fx  ♦ | Mr  (x  - pp) 2 1 0 

(Note:  Equation  8. a is  the  same  as  5. a) 
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The  complexity  of  the  control  equation  increases 
as  we  deviate  from  a constant  forcing  function.  In 
the  main  part  of  the  analysis  we  approximated  the 
forcing  function  as  two  linear  segments. 


F 


Range 


To  calculate  the  potential  as  a function  displace- 
ment for  "approximation  (3)"  would  introduce  too 
much  complexity  into  the  control  circuitry.  The 
only  candidate  control  functions  will  be  approximations 
(1)  and  (2).  (i.e.,  equations  (5. a)  and  (7. a)) 
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(9. a) 

(10. a 


(11. a 
(12. a 


In  terms  of  the  above  parameters,  equations 
(5. a)  and  (7. a)  become. 


-CF< 


x + 


1 

1 


Mr  ‘ m-)2  i 0 


1 

2 


(x 


< 0 


We  now  consider  the  corrections  required  if  other 
forces  act  on  the  recoiling  mass  other  than  the  recoil 
force. 

TfiT, o)  - Kinetic  and  potential  energy  when 

tire  1 ire  weap0n  is  fired 

V - Potential  energy  at  sear  position 

(Kinetic  energy  is  zero) 

W^  - Work  done  by  forces  other  than 

recoil  force  in  going  from  the 
firing  position  to  sear 

From  energy  equation: 


W 


f s 


Vs  - Vfire  * Tf ire 


We  compute  wxs  + v + T until. 


W 


xs 


♦ V + T < T 
— s 


To  determine  the  time  of  fire. 

W - Work  done  on  mass  from  forces  other  than 
the  recoil  force  in  going  from  x to  sear. 
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If  the  recoil  force  is  a constant,  then  the 
control  equation  is  fire  if; 

03-a>  WxS  - Fx  * TMr  <*  - i 0 


The  most  significant  force  would  be  the  weight 
component  when  the  weapon  is  elevated. 

(14. a)  Wxsi  = Mr  g Sin  9 x 


x 


From  equations  (13. a)  and  (14. a) 

(IS. a)  (Mr  g Sin  6 - F)  x ♦ \ Mf  (£  - jjj-) 2 < 0 

For  a 200  lb  weapon  and  90°  elevation  F * 1200  lb. 
The  required  correction  for  the  coefficient  of  x would 

be  ^F  which  is  quite  significant.  Smaller  corrections 

for  friction  and  stripping  of  a round  might  be  included. 
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(16. a 


(17. a 


We  now  investigate  the  correction  due  to  an 
accelerating  reference  frame. 


The  control  equation  is  the  same  as  equation 
(12. a)  except  we  must  include  in  the  term  Wxs  the 

work  done  by  the  inertial  force,  -Mry.  To  see  this 
we  derive  the  energy  equation  for  a moving  reference 
frame . 


/Mr  (x  + y)  dx  = /Fdx 


F = F + F 

cons  neons 


cons 

F 

neons 
AT, 


(X) 


AV 


- Conservative  force  (function  of  x) 

- Non-conservative  force 

- Change  in  kinetic  energy  relative 
to  x-axis 

- Change  in  potential  energy 


From  equation  (16. a)  we  have 


AT 


(X) 


-AV  ♦ / (F „ - M„  y)  dx 

v neons  x ' ' 
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The  term,  ^(Fnc0ns  * Mr  dx»  corresPonds 
to  what  we  called  W , previously.  The  control 
equation  is,  therefore,  fire  if: 


W * T r , 
xs  (x) 


♦ V 1 Vs 


T 


(x) 


Kinetic  energy  relative  to  the  turret. 


If  y is  a constant  and  we  correct  for  elevation, 
and  the  recoil  force  is  a constant,  the  control 
equation  is: 

fire  if, 


(18. a)  (Mr  g Sin  6 - F - Mfy)  x ♦ \ Mr  (x  - 2 < 0 


For  a 200  lb  weapon,  and  1-g  acceleration, 
and  F ■ 1200  lb  the  correction  to  F would  be 

hr. 


Finally,  we  consider  a time  delay  that  is 
constant  from  round  to  round.  Because  of  the 
self -correcting  nature  of  the  control  equation, 
small  variations  in  time  delay  will  not  be 
detrimental . 

T a , - Time  between  firing  and  when  the 

el  ' round  effectively  delivers  the 
impulse  to  the  weapon 

To  determine  the  firing  time  we  must  project 
ahead  the  displacement  and  velocity. 
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(20. a) 


(21 . a)  x 


x , x„  - measured  displacement  and  velocity 
m ’ m r ' 

F - external  forces 
e 

t - time  of  measurement 

m t T , 

m + delay 


m 


dt 


t + T 
m delay 


m 


m delay 


x_  ♦ x„  Tj_i  + / j/  Ctm  ♦ T,a1,„  - T)  dT 

r 


m 


delay 


t 


m 


For  this  correction  we  assume  all  forces  are 
constant  (recoil,  gravitational,  and  inertial). 


(22. a) 


Fg  = F - Mr  g Sin  e - Mr  y 


To  correct  for  ignition  delay  we  add  to  the 
measured  velocity  and  displacement  at  a particular 
time  - a constant. 


(23. a) 


• • 0 

X = X ♦ rr-  I , - 

m Mf  delay 


+ Fe  ^delay)2 


x = xm  + xm  Tdelay  T fT  7 


25 


In  the  computer  we  allow  for  making  or  not 
making  corrections: 


* tie  lay 


I 

e 


I 

P 


1 

0 

1 

0 

1 

0 

1 

0 


If  we  correct  for  time  delay 
Otherwise 

If  we  correct  for  elevation 
Otherwise 

If  we  correct  for  inertial 
forces 

Otherwi se 

If  we  use  a constant  function 
to  compute  spring  potential 

If  we  use  a linear  function 
to  compute  spring  potential 


The  control  numbers  allow  us  to  determine  the 
effect  on  the  system  if  a correction  is  made  or  not. 


(25. a) 


Summary : 

F_  ♦ F. 


1 le  Mr  8 Sin  9 - h Mr  X 


X+I  (X  T ♦ — ^Tde  \ 

xm  *delay  ^Xm  ^delay  Mr  ‘ 


x = x + I 
m 


0 

delay  Tdelay^ 
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Fire  if. 


(F,  ♦ F ) F - F v2 

(25. b,  - S-j  ?.  X . (1  - . ) (S-j-J.  - Fsx, 


+ (Ie  Mr  8 Sin  6 + Ij  Mr  y) 


* IMr  <*  - r-)J  i » 
r 


What  has  been  called  x in  the  previous  analysis 
is  the  center  of  mass  of  the  system  with  respect  to 
the  turret.  If  the  moving  parts  inside  the  weapon 
are  light,  then  we  can  measure  the  receiver  displace- 
ment and  velocity,  otherwise  we  need  to  measure  the 
center  of  mass  of  the  system. 


Qcm  - Displacement  we  desire  the  center 
of  mass  returned  to  after  firing. 


0Cmo  " Velocity  of  the  center  of  mass  we 
desire  the  center  of  mass  returned 
to. 


(26. a) 


- Q 


cm 


O 


- Q 


cm 


O 


For  example,  if  we  start  with  the  slide  and 
receiver  at  sear  with  zero  velocity,  then 


q - Position  of  receiver  at  sear 
1 o 


q - Position  of  slide  at  sear 
2 0 
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Position  of  chamber  at  sear 


(27. a)  Q = (q  W + q W . + q „ W .)  / (W„„„  ♦ 
v ’ xcm0  Min  rec  M2o  si  co  c*'  rec 


W , + W . ) 
si  ch'' 


Q =0 
xcme 


Since  the  chamber  is  always  locked  to  the 
receiver  when  the  weapon  fires,  we  don't  need  to 


rec  T "si  T "ch'  "i  T "si  V'  '"rec 


measure  qc- 

Q 

xcm 

[fW  + w 

1 v rec 

W , + W , ) 
si  ch-7 

• 

Q * 
xm 

[ (W  ♦ W 

W . + W , ) 
si  ch-' 

Wci  ♦ W .)  ,,  ♦ W - qJ/(Ww  + 


From  equations  (27. a)  and  (28. a), 


* q20  ' q10)  "»1  * ‘“CO 


(29. a)  x 


(q  -q  ) ♦ 

i i o 


w ♦ w , ♦ w . 

rec  si  ch 


WS1  q: 


q ♦ 

1 W ♦ W . ♦ W . 
rec  si  ch 
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From  equation  (29. a)  we  see  the  second  term 
gives  the  error  in  measuring  only  the  receiver 
displacement  and  velocity  as  a measure  of  the  center 
of  mass  displacement  and  velocity. 


W + 
rec 


W . + 
si 


W 


ch 


200  lb 


Approximate  values: 


w . = 

si 

32 

lb 

Wch  = 

16 

lb 

q10  = 

0 

q20  = 

9 

in 

qco 

7 

in 

When  the  weapon  fires, 


q 


2 


0 


Error  in  displacement  at  time  of  fire: 


9(32)  +.7C16)  = 2 

200 


For  five  inches  of  recoil  the  error  would  be 
40%  which  is  not  tolerable.  Also,  computer  results 
indicate  that  the  impact  of  the  slide  into  the 
receiver  would  give  a velocity  to  the  receiver 
large  enough  to  cause  firing  of  the  weapon.  Conclusion 
the  receiver  cannot  be  used  as  a measure  of  the  C.G. 
of  the  system. 
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On  page  119  there  are  two  types  of  time  delay; 
ignition  delay  (t^elay^  and  a delay  (t()  because 

of  the  round  impulse  being  distributed  over  time 
rather  than  being  a delta  function.  We  now  show 

how  to  determine  t . 

l 


C (t,tc)  - Chamber  force 

tj-  - Firing  time 

t - Cycle  time 

c 

F - External  forces  acting  on  system 

other  than  the  chamber  force 


C (t,tc) 


Mr  x - F - C (t,tf) 
For  one  cycle: 


(30. a)  0 


I 


f 


(F  - C(t ,tf ) ) dt 
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(31. a 


(32. a 


(33. a 


(34. a 


(35. a 


(36. a 


0 


Mr  x 


(F  - C(t,tf))  dt  dt 


From  equations  (30. a)  and  (31. a) 


C(t,tf))  t dt 


0 


We  compute  the  cycle  time  assuming 
c(t,tf)  = 1 6 (t  - t£)  and  F = constant. 


t 


c 


I 

F 


A nominal  firing  time  is  computed  based  on 
the  same  assumptions. 


tf(n)  " TF 

From  equation  (32. a)  and  assuming  tc  is  a 
constant  given  by  equation  (33. a) 

F t 2 *c 

— - j C (t , tf ) tdt  - 0 

2 


tf(n)  ' tf  " ti 


Equation  (36. a)  represents  the  time  delay  in 
the  system  due  to  the  distributed  force,  c(t,t£). 
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We  now  show  how  to  calculate  tf  from  equation 
(35. a).  t 


C(t,tf) 


tf  t4  ti^l 


(37. a)  C(t,tf)  - ai  (t  - t£)  + b£ 


^ 1 1 ti  ♦ 


1 


From  equations  (35. a)  and  (37. a)  we  have 


F t 2 A.  b. 

C38.a)  l (t.n2.  t.2,  . -A  (t 


2 . t 2) 
i+1  *i  > 


a. 


7 (Vi3  - *i3> 


32 


(39. a) 


(40. a) 


h. 


a . 


r-  <4.i?  - 4?)  + r <4 


i . 1 


4’> 


Ft  2 
c 


a . 


2~ 


- (t.  2 - 

l4  + l 


V> 


From  equations  (33. a),  (34. a),  and  (39. a) 


bi 


a.- 


ti  = 


I 

TV 


^i+l2  * V)  + 31  (4.1 


3 . 


I2 


V>  - IF 


a . 


I1  (4 


2 - t 2i 

i+ 1 4 ] 


C...  - C. 
Note:  a.  = - 1 —■ 

1 4 + l ' 4 


and 


bi  = 4 - 4 4 


See  special  computer  program  for  computing  ti. 

tj  = 10  m sec,  computed  for  chamber  force  with 
I = 160  lb/sec 
F - 1200  lb 

and  using  supplied  C ( t , t £ ) . 
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(reverse  blank) 


CHAMBER  MOMENT  OP  INERTIA 


Moment  of  Inertia  Calculation  for  Chamber. 


(I) 


l22  * (7  mb2  ' 7 ma2)  lb.S^C 

m2 

Density  x Volume 
m = * — — 

g 

g * 32.17  ft/sec2  = 386.04  in/sec2 
Vol  - u (b2  - a2)  L 

m (386.'Oi)  Vo1 


aI 


Vol  r 


4.  36 
~T~ 

3.40 
— T~ 


2.18 

1.70 


1.35 

tt ( 2 . 1 8 2 - 1.702)  1.35 


7.8987 
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0.002885 


(II) 


M, 


1 

2 


,0.282 

'•386.04 


) Vo  1 


I 


accounts  for  lugs 


Wj  = Mj  g = 1 .1137167  lb 


! 2 2 = y~  ( 2 . 1 8 2 - 1 . 7 2 ) = 0.00268651 


lII 


3.40 

~~2 


1 . 26 
~2 


1 . 7 
0.63 


Ln  = 1.09 

Voln  = n ( 1 . 7 2 - 0.632)  1.09  = 8.53721 

Mu  * <TS 57OT>  Vo1  • '>■006236 
W„  = Mn  g = 2.4075  lb 
I77  = ( 1 . 7 2 - 0.632  ) = 0.007773 

ZZji  £ 


(III) 


3.4 


’III 


= 1.7 


al  II  = = 1,026 


L 


III 


7.53 
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Using  Prismoidal's  formula  we  find 

Volt  = L (ri2  + rir2  ♦ r22)  = 2.040818 

Ml  = (3^-3- 4)  Vo  1 1 = 0.0014908 
Wl  = Mi  g = 0.5755  lb 


L2  = 0.5 

Vol 2 = nb  2 2 L 2 = 4.5396014 

m2  = (^g-g-^4)  Vol2  = 0.00331615 
W2  = M?  g = 1.280168  lb 

1 22z  = I M2  b22 
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22 


IV 


1 2 2 •>  ' 1 2 2 


I,,  = \ M2  (1.7)2  - ffl-  M (1.252  . i 


22  ' 7 2 


I = 0.004792  - 0.000228  = 0.004564 

^IV 


1 22  = 1 2 2 j + 1 2 2 j j + I22ni  + X22JV 


I22  = 0.044191 


w = wi  + wn  + wm  + wIV 


W = 16.4828 


0 2 6 2 ) 
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OUTPUT  DATA  (TIME  DEPENDENT) 


1 . 

TIME 

Time  (sec) 

2. 

Q(D,  qi 

Receiver  displacement  relative  to 
turret  (in) 

3. 

Q C 2 ) , q 2 

Slide  displacement  relative  to  the 
receiver  (in) 

4. 

QD(1),  q! 

Receiver  velocity  relative  to  turret 
(in/sec) 

5. 

QD  (2) , q2 

Slide  velocity  relative  to  the 
receiver  (in/sec) 

6. 

QDD(l) , qj 

Receiver  acceleration  relative  to 
the  turret  (in/sec2) 

7. 

QDD(2) , q2 

Slide  acceleration  relative  to  the 
receiver  (in/sec2) 

8. 

pc,  pc 

Cylinder  pressure  (PSI) 

9. 

TC,  Tc 

Cylinder  temperature  (°R) 

10. 

TCD , Tc 

Time  derivative  of  cylinder  temperature 
(°R/ sec) 

11. 

WO,  Wo 

Flow  rate  through  orifice  ( * b-j  -5— ■ ) 

12. 

PSL,  PM 

Barrel  (receiver)  pressure  (PSI) 
(build  up  begins  when  projectile 
passes  port) 

13. 

TBL  * Tbl 

Barrel  (receiver)  temperature  (°R) 

14. 

XMACH,  Mo 

Mach  number  at  throat  (dimensionless) 

15. 

WOS 

Total  mass  flow  (/wQdt)  (lb~?ec2) 

in 

16. 

GAS  IMP 

Total  impulse  applied  to  slide  from 

cylinder  (fpc\.dt)  (lb-sec) 
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17. 

QC.  qc 

Chamber  displacement 

18. 

QCP 

(dimensionless) 

19. 

QCPP 

(in'1 

dq22 

20. 

QCD,  qc 

Velocity  of  the  chamber  relative  to 
receiver  (in/sec) 

21. 

THETC,  ec 

Rotation  of  the  chamber  (rad) 

22. 

THETCP 

(aff)  (rad/in) 

23. 

THEPP 

(rad/ in2) 

dq22 

24. 

THETCD 

Angular  velocity  of  the  chamber  (f^r) 

1 25. 

FR(  1) 

Recoil  Force  (lbs) 

26. 

FR(2) 

Viscous  damping  action  between  receiver 
and  turret  (lbs) 

27. 

FR  (3) 

Forward  buffer  force  of  receiver  (lbs) 

28. 

FR  ( 4 ) 

Rearward  buffer  force  of  receiver  (lbs) 

29. 

FR  (5) 

* FS  ( 1 ) (lbs) 

30. 

FR(6) 

-FS (2)  (lbs) 

31. 

FR  ( 7) 

-FS (3)  (lbs) 

32. 

FR  (8) 

-FS (4)  (lbs) 

33. 

FR  ( 9) 

-FS  (5)  (lbs) 

34. 

FR (10) 

-FS  (6)  (lbs) 

35. 

FR  (11) 

-FS (7)  (lbs) 
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36. 

F R ( 1 2 ) 

- FC  ( 1 ) (lbs) 

37. 

FR ( 1 3) 

- FC  ( 2)  (lbs) 

38. 

FR(14) 

- FC ( 3)  (lbs) 

39. 

FR ( 1 5) 

Piece-wise  constant  dissipative 
(coulomb  friction)  (lbs) 

force 

40. 

FS  ( 1 ) 

Drive  spring  on  slide  (lbs) 

41. 

FS  (2 ) 

Viscous  damping  between  receiver 
slide  (lbs) 

and 

42. 

FS  (3) 

Cylinder  pressure  force  (lbs) 

43. 

FS  (4) 

Slide  rear  buffer  force  (lbs) 

44. 

FS  (5) 

Sear  spring  force  (lbs) 

45. 

FS  (6) 

Slide  forward  buffer  force  (lbs) 

46. 

FS  (7 ) 

Piece-wise  constant  dissipative 
(lbs) 

force 

47. 

FC  Cl) 

Chamber  rear  buffer  force  (lbs) 

48. 

FC  ( 2 ) 

Viscous  damping  between  receiver 
chamber  (lbs) 

and 

49. 

FC(3) 

Piece-wise  constant  dissipative 
(lbs) 

force 

50. 

FC  (4) 

Propellant  force  (lbs) 

51. 

VP0TQ  (1 ) 

Weight  force  (f^— ) (lbs) 

52. 

VP0TQ  (2) 

Weight  force  (|^-)  (lbs) 

53. 

ENERGY 

Parameter  computed  to  determine 
(in-lbs) 

firing 

54. 

FXR 

Sum  of  forces  acting  on  receiver 

(not 

weight  or  constraint  forces) (lbs) 
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55. 

FXS 

Sum  of  forces  acting  on  slide  (not 
weight  or  constraint  forces)  (lbs) 

56. 

FXC 

Sum  of  forces  acting  on  chamber  (not 
weight  or  constraint  forces)  (lbs) 

57. 

XMC ( 1 ) 

Sum  of  moments  acting  on  chamber 
( in- lbs) 

58. 

TXT 

x-coordinate  of  turret  (center  of 
rotation)  in  inertial  frame  (m) 

59. 

TXTP 

Velocity  (x-coordinate)  in  inertial 
frame  (m/sec) 

60. 

TXTPP 

Acceleration  (x-coordinate)  in  inertial 
frame  (m/sec2) 

61. 

TYT 

y-coordinate  of  turret  (center  of 
rotation)  in  inertial  frame  (m) 

62. 

TYTP 

Velocity  (y-coordinate)  in  inertial 
frame  (m/sec) 

63. 

TYTPP 

Acceleration  (y-coordinate)  in  inertial 
frame  (m/sec2) 

64. 

TZT 

z-coordinate  of  turret  (center  of 
rotation)  in  inertial  frame  (m) 

65. 

TZTP 

Velocity  (z-coordinate)  inertial  frame 
(m/sec) 

66. 

TZTPP 

Acceleration  (z-coordinate)  inertial 
frame  (m/sec2) 

67. 

PSI 

Euler  angle,  T,  (rad) 

68. 

PSIP 

f (rad/sec) 

69. 

PSIPP 

f (rad/sec2) 

70. 

THETA 

Elevation,  0,  (rad) 

71. 

THETD 

0 (rad/sec) 

72. 

THEPP 

0 (rad/sec2) 
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OUTPUT  DATA  (TIME  INDEPENDENT) 


1 . 

THEOC,  9 

oc 

2. 

YOCI,  Y 

’ OCl 

3. 

AC1 , Aci 

4. 

XCI,  xci 

S. 

QCI,  qci 

6 . 

Y°CF,  Yocf 

7. 

THEOF,  0 

oc 

8. 

ACF,  a r 
ct 

9. 

XCF,  xcf 

10. 

VCF,  ycf 

11 . 

QCF,  qcf 

12. 

Y°QF’  Yoqf 

13. 

theoq,  eoq 

14. 

AQF.  aqf 

15. 

UQF,  uqf 

Angle  associated  with  kinematic 
constraints  (rad) 

Distance  associated  with  kinematic 
constraints  (in) 


Major  axis  of  ellipse  for  kinematic 
constraint  (in) 

Distance  associated  with  kinematic 
constraints  (in) 

Distance  associated  with  kinematic 
constraint 

Distance  associated  with  kinematic 
constraint  (in) 

Angle  associated  with  kinematic 
constraint  (rad) 


Distance  associated  with  kinematic 
constraint  (in) 

Distance  associated  with  kinematic 
constraint  (in) 

Distance  associated  with  kinematic 
constraint  (in) 

Distance  associated  with  kinematic 
constraint  (in) 

Distance  associated  with  kinematic 
constraint  (in) 


Angle  associated  with  kinematic 
constraint  (rad) 

Distance  associated  with  kinematic 
constraint  (in) 

Distance  associated  with  kinematic 
constraint  (in) 


45 


16. 

QC0'  ^co 

Distance  associated  with  kinematic 
constraint  (in) 

17. 

RB,  Rb 

in  2 

Gas  constant  ( ) 

sec2  - °R 

18. 

GAMMAB,  Yb 

Ratio  of  specific  heats  (dimensionless) 

19. 

RH0A,  pa 

d 

Density  of  air  (lbs-sec2/in4) 

20. 

CRRAT 

Critical  pressure  ratio  (dimensionless) 

21. 

XMASS(I) 

Mass  of  receiver,  slide  and  chamber 
(lbs-sec2/in) 

22. 

QREF ' Rref 

Position  of  slide  in  equilibrium 
position  (in) 

23. 

RIMPUL 

Chamber  impulse  (lbs-sec)  (/c(t)dt) 

24. 

TOTMAS 

Total  mass  of  receiver,  slide  and 
chamber 

25. 

VFr 

FRECO,  s2- 

Average  recoil  force  (lbs) 

26. 

RECOLD,  L 

Nominal  recoil  distance  (in) 

27. 

YDDOT,  y 

Acceleration  of  turret  to  correct  for 

firing  time  (in/sec2) 
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FLIGHT  PATH  CALCULATIONS 


The  flight  path  is  represented  by  a series  of 
linear,  circular,  and  spiral  segments.  It  uses  the 
same  flight  path  representation  as  contained  in 
"The  FUE  Air  Defense  Evaluation,"  Fred  L.  Bunn, 

Vol  3 of  3.  Excellent  documentation  is  contained 
in  Appendix  A of  the  above  report  of  the  two  routines 
required  to  represent  the  flight  path.  Slight 
modifications  were  made  to  compute  the  acceleration 
along  the  X0Y0Z0  coordinates.  Also,  the  first  segment 
should  be  a line  segment  and  no  two  line  segments 
should  be  used  consecutively,  which  would  imply  an 
instantaneous  change  in  direction  of  the  velocity, 
which  is  impossible. 

What  is  the  XYZ  axis  in  the  "Fue  Air  Defense 
Evaluation"  corresponds  to  Z0X0Y0  in  this  analysis. 

The  vehicle  is  started  at  the  origin  in  the  inertial 
reference  frame.  Angular  rates  of  the  turret  are 
obtained  by  assuming  the  weapon  is  always  pointed 
at  some  fixed  point  in  space.  By  changing  the  point 
in  space,  we  obtain  different  angular  rates  and 
angular  accelerations. 


yo 


47 


6Q  - Heading  angle  of  aircraft  at  t = o 
6 ^ - Heading  angle  at  a later  time 


Tt  - Vector  in  direction  of  target 


(x(t),  y(t),  z(t) ) - Location  of  vehicle 


(xt,  yt,  zt)  - Location  of  target 


- I(x(t)  - xt)  x0  ♦ Cy (t)  - yt)  y0  ♦ (z(t)  - Zt)  z0] 


xt)2  ♦ (y (t)  - yt)2  «■  ( z ( t ) - z t ) 2 


8 


See  main  analysis  for  rotational  transformation. 


From  geometry  we  have 

(3.b)  Tt  = Cos  y Cos  e x0  + Cos  y Sin  0 y0  - Sin  y 


(4.b) 


From  equation  (3.b) 


tan  8 = 


y(t) 

xitr 


- X. 


(5.b) 


0 = tan  1 ( 


y(t)  - y. 


xTtT 


- x. 


(Use  ATAN2  in  computer) 


tan  y = 


(z(t)  - zt)  Sin  0 
(y(t)  - yt) 


From  equation  (4.b) 


(6-b)  ® = (x('0S-  tit*)  (x(t)  ' xt)  - (y (t)  - 


From  equation  ( 5 . b ) 


Cos  2 y 


(7>b)  y ‘ 'fy(t) *yt)2  Hz(t)  Sin  0 + (z  (t)  - z t ) 

(y(t)  - yt)  + ( z ( t ) - z t ) 


yt)  x(t)J 

e Cos  0) 
y (t)  Sin  ej 
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(8.b)  0 * [{-0  Cos  o Sin  0 (y(t)(x(t)  - xt) 

- Cy(t)  - y t)  x(t))  + Cos2  0 (y(t) (x(t)  - xt) 

- (y (t)  - yt)  X (t ) ) } (x(t)  - xt)2 

- 2 Cos2  0 (y (t) (x (t)  - x t ) - (y  (t ) - yt)  x(t)) 
x(t)(x(t)  - xt)J/(x(t)  - xt)4 

(9 . b)  y = -[{-y  Cos  y Sin  y ( ( z ( t ) Sin  0 ♦ (z(t)  - zt) 

0 Cos  0)  (y (t)  - yt)  + C z ( t ) - zt)  y(t)  Sin  0) 

♦ Cos2  y [y(t)(z(t)  Sin  0 + (z(t)  - zt)  0 Cos  0) 

♦ (y(t)  - yt)(z(t)  Sin  0 ♦ z(t)  0 Cos  0 + z(t)  0 
Cos  0 - 02  Sin  0 (z(t)  - z^)  ♦ (z  (t ) - zt)0  Cos  0) 

♦ z(t)  y(t)  Sin  o ♦ (z(t)  - zt)  0 y(t)  Cos  0 

♦ ( z ( t ) - z t ) y (t)  Sin  0J)  (y(t)  - yt)2 

- 2y(t)  Cos2  y (y(t)  - yt)  [(z(t)  Sin  0 + (z(t)-zt) 
0 Cos  0)  (y ( t ) - yt)  ♦ (z(t)  - z t ) y(t)  Sin  0j] 

/ (y(t)  - yt) 4 
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Side  View 


f 


X0  , 

y° , 

Zo 

Inertial 

reference  frame 

*1  , 

yi » 

Zi 

Origin  at 
(parallel 

center  of  rotation 
with  x0,  yQ  , Zo) 

*2  » 

y2 » 

Z 2 

Origin  at 
(rotated 

center  of  rotation 
about  z,  through  angle  0) 

i 


51 


Top  View 


Xf,  y^,  Z£  - Fixed  body  axis  (rotated  about  y2 , 
through  angle  V) 


The  center  of  rotation  is  allowed  to  move  in  space; 
position  is  given  by: 

rcg  = x(t)  x0  ♦ y(t)  y0  ♦ z(t)  z0 

Also,  the  two  "Euler  Angles,"  6 and  v , are  described 
as  functions  of  time  (the  first  derivatives  represent 
slewing  rates) . 
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We  now  derive  transformation  relationship  between 
coordinate  systems: 


x 2 = Cos  e xi  + Sin  9 yi 
y2  = -Sin  0 x i + Cos  0 yi 
z2  = zi 

where  x2  is  a unit  vector  along  x2  axis. 


Let 


(column  matrix) 


Then  equation  (1)  can  be  written  as: 


From  equations  (2)  and  (4) 


Cos  y o 

-Sin  4* 

Cos 

0 S in  0 

0 

(5) 

XI 

II 

0 1 

0 

-Sin 

0 Cos  0 

0 

Sin  t 0 

Cos  4* 

0 

0 

1 

* 

Cos  T Cos 

e Cos 

t Sin  e 

-Sin  V 

(6) 

Xf  = 

-Sin  0 

Cos  0 

0 

X, 

Sin  ¥ Cos 

0 Sin 

V Cos  0 

Cos  t 

We  now  determine  the  "kinematics"  of  the  receiver. 
We  first  determine  the  angular  velocity  and,  secondly, 
we  determine  the  C.G.  relative  to  the  inertial  frame 
of  reference.  (Angular  velocity  is  determined  relative 
to  the  fixed  body  axis,  x^,  y^ , z^) 

The  angular  velocity  of  the  receiver  is  given  by, 

(7)  Wm  = 0 zi  + i y2 

X£  is  of  the  form, 

(8)  Xf  = A Xi 

and  from  equation  (8), 

(9)  Xj  = A* 1 Xf  * AtXf 
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From  equations  (6)  and  (9)  we  have 


(10) 


X,  = 


Cos  y Cos  6 -Sin  0 

Cos  y Sin  6 Cos  0 

-Sin  V 0 


Sin  t Cos  0 
Sin  y Sin  0 
Cos  y 


(ID 


From  equation  (1)  it  can  be  seen  that. 


zi  = -Sin  y xr  + Cos  y z. 


From  equation  (4)  we  have, 


(12) 


Cos  y 0 Sin  y 

0 10 

-Sin  y 0 Cos  y 

From  equation  (12) , 


(13) 


Yz  * yf 


From  equations  (7),  (11),  and  (12)  we  have, 


(14) 


= -e  Sin  y xf  + y yf  + 0 Cos  y zf 
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* 


(x'r»  y\,  z'r)  - Initial  position  relative  to 

xf»  Yf,  Zf  coordinate  system. 

xr,  yr,  zr  - Fixed  body  axis  of  the  receiver. 

( xr  is  parallel  with  translation  axis) 

q l * Travel  along  the  translation  axis  of  the  C.G. 
of  the  receiver. 


5 


The  C.G.  of  the  receiver  relative  to  the 
inertial  frame  of  reference  is. 


(15) 


rmr  = x(0  x°  + y(t)  yo  ♦ z(t)  Zo  ♦ (q  1 ♦ x'r)  x f 


+ y'r  yf  + z r zf 

From  equations  (6),  (15),  and  the  fact  that, 


(16) 


xi 

yi 

Zi 


Xo 

yo 


(17)  rm  = {x(t)  + Ui  * *'-)  Cos  f Cos  9 - y'  Sin  0 ♦ 


z'r  Sin  ¥ Cos  0 } Xo 


+ (y(t)  ♦ (q  i ♦ x 'r ) Cos  ¥ Sin  0 ♦ y'  Cos  0 ♦ 
z'r  Sir,  ¥ Sin  0 } y0 


+ {z(t)  * (qi  ♦ x'r)  Sin  t + z'r  Cos  v}  Zo 
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Also,  from  the  fact  that, 


(18) 


X 


r 


we  can  rewrite  equation  (14)  as, 


(19) 


17  = -6  Sin  V x + ¥ y +0  Cos  f z 

m_  r 7 r r 


The  computer  operates  on  equations  (17)  and  (19) 
to  generate  equations  of  motion;  see  part  A. 

Since  Wm  is  independent  of  the  degrees  of  freedom, 
r 

the  equations  of  motion  will  be  independent  of  Wm  . 

Therefore,  orientation  of  xf , y , zr  is  not  important 

for  determining  moments  of  inertia,  hence  we  do  not 
have  to  calculate  them. 


Equation  (17)  gives  the  C.G.  of  the  receiver 
relative  to  the  inertial  reference  frame.  Equation 
(19)  gives  the  angular  velocity  of  the  receiver 
relative  to  the  xr , y , zr  coordinates.  The 

coordinate  qi  must  be  determined  yet. 


We  consider  forces  applied  to  the  receiver  along 
the  translational  axis  as, 


(20) 


x 


r 
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The  point  of  application  is  given  by, 


(21) 


r ' 


r 


From  equations  (15)  and  (18)  we  have, 

(22)  rm  = x(t)  x0  + y(t)  y0  + z(t)  z0 
r 


♦ (m  + x'r) 


(23) 


The  generalized  force  equation  is. 
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3rr 

if  1 * 1 1^7  = 0 

if  i = 1 from  equation  (22), 


(24) 


3q ! 


From  equations  (20),  (23),  and  (24) 


(25) 


| 


F includes  all  forces  on  receiver  (except 
xr 

gravitational  and  constraint  forces)  along  translation 
axis . 
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Dynamics  of  Internal  Parts 


Receiver 


Chamber 


ys  t. 


Slide 


xc  ^c  zc 


<xv  yv  zV 


(X'c»  y'c*  z'c) 


mass  of  slide 
mass  of  chamber 

coordinate  system  fixed  to  the 
chamber  (coincides  with  principal 
axis) 

slide  position  at  t = 0 relative 
to  xr  yr  z r coordinate  system 

chamber  position  at  t * 0 relative 
to  xr  yr  zr  coordinate  system 


N°te:  xs  = xr 


The  angular  velocity  of  the  slide  is  given  by: 


* -e  Sin  y x_  ♦ y y_  ♦ 0 Cos  y z 


The  C.G.  of  the  slide  relative  to  the  inertial 
frame  of  reference  is, 

<27>  Fm  ‘ Fm  " (x's  * ‘f  J xr  * yr  * z's  zr 

s r 

From  equations  (15)  and  (18)  we  can  rewrite 
equation  (27)  as: 

(28)  F = x ( t ) x0  + y(t)  y0  + z(t)  Zo  + (q j + q2  + x'<- 

in  ** 

s 

+ X'r)  xf  + (y's  + y-r)  yf  + (z'r  + z's)  zf 

From  equations  (6)  and  (16)  we  can  rewrite 
equation  (28)  as: 

(29)  rm  = [x(t)  + (qx  + q2  + x'x  + x'r)  Cos  y Cos  0 

ms 

- (y ' + y'r)  Sin  0 + (Z'r  + Z's)  Sin  y Cos  0]  xe 

+ t y ( t ) ♦ (q  j + q2  + X's  + X'r)  Cos  y Sin  e 

+ (y-  + y'r)  Cos  o + (z'r  + z-s)  Sin  y Sin  0]  yc 

+ [ z ( t ) - (q  j + q2  + x's  + x'r)  Sin  y 

+ (z'r  * *'*>  Cos  Y J Zo 
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The  computer  operates  on  equations  (26)  and  (29) 
to  obtain  expressions  for  equations  of  motion. 


We  shall  be  interested  in  the  generalized  forces 
from  the  force  applied  to  the  slide  along  the 
translational  axis  at  the  point  r . 


r 


From  equation  (28)  and  the  fact  that  xf  * x 

I s 
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Also, 


(31) 


Q? 


3q2 


F 

x 


s 


From  equations  (30)  and  (31) , the  generalized 
forces  are  independent  of  the  point  of  application 
on  the  slide. 


We  now  obtain  the  kinematics  of  the  chamber. 

The  rotation  of  the  chamber  is  a function  of  q2 
which  we  call  0 (q2).  The  translation  of  the  chamber 

has  a certain  amount  of  motion  relative  to  the  slide 
when  the  weapon  is  unlocked,  but  we  will  not  add 
another  degree  of  freedom  to  account  for  this  relative 
motion . 

The  location  of  the  C.G.  of  the  chamber  is  given  by 


(32) 


m 


= r 


m 


(x'. 


qc) 


+ y y + 
/ c 7 r 


c zr 


(qc  is  a function  of  q2) 

Equation  (32)  is  of  the  same  form  as  equation 
(27)  with  q 2 replaced  by  qc  and  x'g,  y'g,  z' 

replaced  by  x 'c , y'  , z'  respectively.  Therefore, 

from  equation  (29)  we  can  rewrite  equation  (32)  as: 
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(33) 


[x(t)  ♦ (q  i ♦ qc  ♦ x 'c  + x"r)  Cos  y Cos  e 


- (y'c  ♦ y\)  Sin  0 ♦ (z'r  ♦ z'c) 


Sin  y Cos  ej  xQ 


+ Iy(t)  ♦ (q i ♦ qc  * x'c  ♦ x'r)  Cos  y Sin  e 

* (y'c  + y'r)  Cos  0 + (z'r  + *-c) 

Sin  y Sin  0]  y0 


+ £z(t)  - (q  i * qc  ♦ x'c  ♦ x'r)  Sin  y 


♦ (z^  + *'c)  Cos  y]  z0 


For  forces  acting  along  the  translation  axis  on 
the  chamber,  the  generalized  forces  are  given  by, 


(34) 


Qi 


(35) 


Q2 


3q2 


66 


The  angular  velocity  of  the  chamber  is  given  by, 

(36)  \ - v * e tq2>  *c 

c r 

• From  equation  (19)  we  can  rewrite  equation  (36) 

as , 

(37)  » -0  Sin  xr  + $ yr  + e Cos  v zr 
mc 

+ ®c  (q2)  xc 


(38)  xr  - xc 

yr  - Cos  ec  yc  ♦ Sin  ec  zc 
zr  - -Sin  ec  yc  ♦ Cos  ec  zQ 
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From  equations  (37)  and  (38)  we  have, 


(39)  =■  [0C  (q2)  - 6 Sin  vj  xc 

c 

♦ Cos  ec  - e Cos  Y Sin  0cJ  yc 

♦ [Y  Sin  0C  + 0 Cos  Y Cos  0CJ  zQ 

5c  - Hi  ** 

The  generalized  forces  due  to  a moment  M applied 
to  the  chamber  is  given  by:  c 

<41>  Qi  - Mc  - 0 

t42)  Qa  - Mc  ||£ 

To  derive  equations  (41)  and  (42)  we  have  applied 
a result  that  is  available  in  Part  A,  page 

The  gravitational  potential  is  given  by: 

V ' Mr*  (Fmr>  * Mcg  (Fm  ) 

/-component  c /-component 

* Ms*  (Fm  > 

y-component 
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With  no  vehicle  motion: 


xt  yt 


Kinematic 
Fm_  " 


z - coordinate  system  in  direction  of 
translation  axis 

relationships  are: 

(qi  ♦ *-r)  xt  ♦ y'r  xt  + z'r  zt 


0 


A A 

(qi  ♦ qc  ♦ x'r  + x'c)  xt  + (y ' T ♦ yc)  yt 


♦ 


z 


t 


W * 0 (q2)  x„ 

m c c 

c 

V - Mr  g Sin  0 (qi  + x'r) 

Mg  g Sin  0 Cq l ♦ q2  ♦ x'T  * x's) 
* Mc  g Sin  6 (q i ♦ qc  ♦ x'r  ♦ x'Q) 


(Generalized  forces  are  the  same  as  the  general 
analysis . ) 
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Locking  and  Unlocking  Kinematics 


To  avoid  discontinuous  changes  in  velocities, 
we  approximate  all  "cam  constraint  functions"  by 
functions  with  continuous  derivatives  (see  receiver 
and  slide  drawings  for  nominal  values). 

The  rotation  of  the  chamber  is  given  by  the 
following  as  a function  of  the  displacement  of  the 
chamber  position  relative  to  the  receiver. 


qc^L  - value  of  q£  when  locking  begins 
qc£k  ' value  of  qc  when  weapon  is  locked 
qci  - location  of  end  of  linear  segment 
qc£  - location  of  start  of  linear  segment 
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We  choose  ellipses  at  the  end  of  the  linear 
segments  to  obtain  continuous  derivatives. 


(1) 


The  equation  of  the 
y = x 

The  general  equation 


line  above  is  given  by 
tan  0 + y0 
for  an  ellipses  is. 


(2) 


x 2 y 2 

— ♦ I—  = 1 


From  equations  (1)  and  (2)  we  have, 


(3) 


dy  _ b2x 

a2y 


tan  0 
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5! 


From  equations  (1),  (2),  and  (3)  the  following 
equations  must  be  satisfied. 


(4) 


x2  + (x  tan  e + y0)2  _ x 
a2  b2 


(S) 


-b2x 


a2  (x  tan  0 + y0) 


= tan  e 


We  now  solve  for  "a"  and  "x"  and  allow  "b"  to 


vary. 


From  equation  (5)  we  have. 


(6) 


h 2 

(—  + tan2  e)  x * -y„  tan  0 
a2 


(7) 


-y0  tan  e 
b2 

- + tan  e 


Let 


(8) 


s - I 
a2 


Substituting  equations  (7)  and  (8)  into  (4) 
we  get, 


(9) 


s 


( -y°  tan  0 + 1_ 

sb2  ♦ tan2  e b2 


C ta»2  L + yo)2 

sb2  + tan2  e 


1 
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* 


(10) 


(11) 


(12) 


(13) 


1 2 
S(-y0  tan  0) 2 + — [-yc  tan2  0 + y0  (sb2  + tan2  0)] 
b 2 

= (sb2  + tan2  0) 2 

S(-y0  tan  0)2  + — S^-  ^ — = s2b4  + 2sb2  tan2  0 + tan4  0 

b2 

S2(b4  - y02b2)  + S[2b2  tan2  0 - (y0  tan  0)2] 

+ tan4  0 * 0 

S = ~f2b2  tan2  e - (yQ  tan  0)21 
2 (b4  - y02b2) 


(14) 


± [2b2  tan2  0 (yQ  tan  0)2]2  - 4(b4  yc2b2)  tan4  0 

2 (b4  - y02b2) 

S = tan2  0 -tan2  e 

y02  - b2  ’ b2 


To  make  "S"  positive  we  use  the  minus  sign  for  the 
discriminat  in  equation  (13),  thus 


(14a) 


S = tan2  Q 
Xo2  - b2 
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From  equation  (8)  , 


(14b) 


Using  equations  (14a)  and  (14b)  the  computation 
of  "a"  and  "b"  is  arbitrary. 


Applying  the  results  from  pages 


Nominal  values: 

« - * 
c F 

«CFL  * 0 
^CIL  = 1-02  in 
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(15) 


>c  = tan'1  ( 


c ■) ' i + n 
^cfL  ^ciL  ^ 


(16) 


y°ci  ~ bci  + eci 


Nominal  values: 


Eci  = 0.05  rad  (chosen  arbitrarily  small  to 
represent  sharp  corners) 

bci  = 0*0^92  rad 


From  equations  (14a)  and  (14b)  with, 


(17) 


yo  = y°ci 

0 = 0o 

c 
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We  compute  "S  ."  from  equation  (14a)  and  "a 
from  (14b).  C1  C1 


(18) 


ci 


From  equation  (7)  we  have, 


(19) 


x . 


Cl 


~y°ci  tan  0 o c 
b .2 

— C 1 + tan2  e o 

a 2 c 


(20) 


lci 


^ciL  + 


X 


ci 


We  will  now  look  at  the  path  prior  to  being 
fully  locked. 
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(21) 


(22) 


(23) 


(24) 


yocf  “ ^b 


c i 

cf  tan  0, 


From  equations  (14a)  and  (14b)  with, 


b = bcf 
Xo  = yocf 


= - (b 


Cl 


cf  tan  e0 


e = ec  , = 0 


cf  ' Doc  ' 7 


We  compute  "S  " from  equation  (14a)  and  "a  " 
from  (14b)  ct  cf 


acf  = / Scf 


Nominal  values: 
bc£  = 0.10  in 

From  equation  (7)  we  have, 


cf 


y°cf  tan  0 ° c f 


bcf2 


cf 


+ tan2  0O  c 

2 Cf 
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UFOR  = (q2  ' qc^  FOR  ‘ value  0f  Cq2  ' qc}  When 

slide  is  fully  forward. 

UpL  = (q  - qc) pl  ' value  of  (q^  - qc)  when 

chamber  is  fully  loaded. 

U I l - (q  - qc)IL  * value  of  (q  - qc)  when 

chamber  is  starting  to  lock. 
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When  the  chamber  is  fully  unlocked  there  is  some 
play  between  the  slide  and  chamber.  This  is  not 
considered  in  the  model  since  adding  it  to  the  model 
would  not  contribute  much  to  understanding  the  response 
of  the  system  and  any  results  obtained  would  be 
questionable  because  there  would  be  an  impact  of  a 
metal  surface  (cam  face)  with  the  roller  which  is 
difficult  to  handle  due  to  the  complexity  and  lack 
of  understanding  in  that  area.  When  the  chamber  is 
unlocked  from  the  receiver  it  becomes  locked  to  the 
slide  and  the  chamber  and  slide  act  as  a single  rigid 
body. 
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(q2  • qc^q£  " value  of  Cq2  - qc)  at  the 
start  of  linear  segment 


angle  "L"  makes  with  the  (q  - q„)  axis. 


(25) 


The  values  of  e^£  and  are  arbitrarily  chosen 
to  provide  a sharp  corner. 


From  equations  (14a)  and  (14b)  with, 


Yo  = Y°qf  - - (bqf  + eqf) 


8 ' 9-«  ' tan'‘  * * 


b = b 


qf 


we  can  compute  "sqf"  and  "a^". 


Nominal  values: 


bqf  = 

0.0344  rad 

Cq£  " 

0.09 

rad 

UFOR  = (q2 

qc)FOR  = °*° 

UFL  = 

(q  2 - 

qc)FL  = *1*25  in 

UIL  " 

(q2  - 

qc ) I L = -2-50  in 
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(26)  Uqf  - (q?  - qc)qf  = UfL  - xqf 

From  equations  (19)  and  (25)  we  have, 

Y f tan  0 „ 

(26a)  Uf  = UCI ^ ^ — 

qf  FL  b 2 

— + tan2  9 
a2  °q 

qf 

We  have  defined  the  two  constraint  functions 
6C  * g (q2  - qc)  and  0£  = f (qc)  which  determines 
qc  as  a function  of  q^.  Since  these  two  constraint 
functions  are  not  simple  functions  we  cannot  expect 
qc  = h (q2)  to  be  simple.  We  first  obtain  a general 

understanding  of  h (q  ) from  f (q  ) and  g (q  - q 1 , 

2 c 2 c 

then  we  show  how  to  compute  h (q  ) along  with  the 

2 

first  and  second  derivatives  which  are  required  to 
analyze  the  dynamics  of  the  system. 


82 


We  use  the  Newton-Raphson  method  to  solve  for 

h(qc)- 

(27) 

ec  = f(qc)  - g(q2  * qc) 

(28) 

qc  » hCq  ) 

To  find  h (q  ),  for  every  value 

2 

find  qc  such  that. 

of  q we  must 

2 

(29) 

G(qc)  - g(q2  - qc)  - f(qc)  = 0 

G'(qc)  = -g  (q2  - qc)  - f (qc) 

Newton-Raphson  algorithm  is, 

(30) 

G(qc(l-l)) 

qc»)  Oc(l-l)  G (q^.^, 

(30a) 

, e'q,  - qc(i-i)) 

' ftqc(i-l)) 

Gc(i,  qc(i-D  g'(q2  - qc(i.1), 

* £’(qc(l-l)> 

(31) 

qCIL  * qCFL 

qc(o)  * 

2 
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(32) 


(33) 


(34) 


(35) 


We  iterate  equation  (30a)  until, 


[g  - £]  < ec 


Nominal  value: 

e = 0.00001  rad 
c 


We  now  obtain  the  first  and  second  derivatives 

of  h(q  ) . 

2 


dqc  dqc  dec  - , 

3q-  = 10-  3q"~  = h (V 

n2  C ^2 

Also  from  equation  (27), 


ec  - f(qc) 


de. 


If  (qJJ 


3^  - 3^ 


f (qc) 


dqc 


f (qc) 
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(36) 


6C  • g(q2  - qc) 


let  x - q - q 
2 *- 


1 


ST  l«C*JJ 

2 


37  igCx)J  3q~ 

2 


(37) 


d0r  - ^c, 

3q“  " g Cx)  U - 3q“] 


From  equations  (33),  (35),  and  (37), 


(38) 


g (q„  - qc)  - 

h (q  ) - 2 [1  - h (q  )] 

2 £ Cqc)  2 


(39) 


h (q  ) 


g (q„  - qc) 
2 L_ 


f (q,.)  ♦ g (q  - qr) 


* g (q  - qc) 

37—  [h(q)J*h  (q  ) * 37—  [— 

Jq2  2 2 dq2  f (qc)  + g (q?  - qc) 


J 


h (q  ) 
2 


[f  (qc)  + g (x)]  jjjj  [g  (x)J  - g (x)  3^J  [f  (qc)+g  (X)  ] 
[f'(qc)  * g (x)]2 
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3q7  [f  (qc)J  " 3^  *f  (qc)J  3q7  = £ h (q2) 

3q7  [g'(x)]  = 37  3^7  = 8"(x)  (1  - h'(q2)) 

(40)  h"  = S U.  .-— h_J.  . g If b ♦ g"  (1  - h )] 

(f  + 8 ) (f  + g )2 


We  are  now  going  to  obtain  expressions  for 
f'(qc)>  f"(qc),  g'(q2  - qc).  and  g"(q2  - qc)  ■ 

For  qc  L qcpL  : 


(41) 


f(qc) 


f'(qc) 


undefined 


f"Cqc) 


For  qCF  < qc  < qCFL  : 


f(qc)  - 5C  * XCF 


‘CP’  'VP2  XV' 

b CF 


qc  “ yCF  + (qCFL  ’ hcF^ 
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l2CF  " 77^"  (qc  ' (qCFL  ' bCF))2jl5 
b CF 


f'(qc) 


(qc  ' qCFL  * bCF) 
f(qC)  • 6c 


f-Cqc)  = I “7  + C- 

fCqc)  ‘ 6C  b2CF 


f'(q.) 


For  qCI  < qc  1 qCF  : 


f(qc)  * Cqc  - qciL^  tan  e0c  + Eci 


(41b)  f ' (qc)  - tan  eoC 


f"(qc)  - 0 
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(42) 


(42a) 


For  U > (q2  - qc)pL  : 


g(u)  = - J 


g'(u)  = 0 
g"(u)  = 0 

For  (q2  - qc)qf  ( u 1 (^  ~ ^FL  : 


g(u)  = bq£  ' l ' ( b 'qf  * ( 


Si 

lqf 


g'(u)  - (g^) 


bqf 2 (U  ~ UFL) 


q£  g(u)  - bq£  +5 


- b‘ 


g"(u)  * 


rL 


[- 


g(u)  - bqf  ♦ £ fl2q£ 
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For  u 1 (q2  ~ q^qf  : 


gCu)  = (u  - ulL)  tan  0 
(42b)  g'(u)  = tan  6 

oq 

g " (u)  - 0 

From  equations  (36)  and  (37), 

ec  - g(q2  - qc) 

de 

(43)  = g'(q2  * qc)  [1  - h. ' (q 2 ) ] 

d2ec 

We  also  need  to  know  

dq  2 2 


(44) 


g " (u)  [1  - h' (q2) J 2 + g ' (u)  h 


(q  2 ) 
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RS(qj) 


p3REF  ' FzREF 
q3REF  * q2REF 


q 3REF^  + FsREF 


Nominal  Values: 

CqiREF»  FlREF^  = in»  1266  lbs) 

(q2REF*  F2REF)  = f2*85  in»  1224  lbs) 

(q3REF>  F 3REF^  ~ C-**7  in,  1136  lbs) 

Calculation  of  recoil  force  from  the  force  deflection 
curve . 


F 


Sear 

Position 
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A l + A 2 * j Mr  (2^)? 

. . 1 I2 

Ai  + A2  = f 


If  A! 


I2 

*5 


then 


qR 


< q2REF 


otherwise , 


qR  - q2REFF 


A i 


q2REF  IF2REF 


Cl  q2REF, 

2 J 


If  Ax 


> 1 12 

* *5 


then 


,CiqR 

qR 


(F2 


REF 


Ci  q2REp)) 


I2 

MR 


,F  2REF  ' Cl  q2REF 
qR  = • C Cl 
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If  Ai  if  RT  » then 


Ai  + (qR  - q2REF)  i F 2REF  + j c2  (qR  - ^2REF)] 


1 i2 
8 M 


R 


p2REF  " I r 2REF 


qR  “ q2REF  + 


[F 


+ F2RPp  - 2C2  (Aj 


2.  Receiver  Friction: 

Assume  that  all  surfaces  moving  relative  to  each 
other  are  well  lubricated.  We  can  assume  under  these 
conditions  that  dissipative  forces  acting  on  the 
receiver  are  dependent  on  only  the  velocity  relative 
to  the  turret. 


FR2  = * CRz  qi 
where 


CR  - Frictional  coefficient 
2 (Nominal  Value  = 0) 

3.  The  forward  and  rearward  buffers  are  assumed  to 
be  linear  springs,  and  the  energy  absorbed  during 
compression  and  decompression  is  assumed  to  be  due 
to  "hysteresis."  The  energy  loss  is  represented  to 
be  some  percentage  of  the  total  potential  (when  mass 
has  come  to  rest)  of  the  spring. 
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9 


F 

rPLD 

ELOSS 
t loss 

q 

o 

qi 

q' 


Preload  force 
Energy  loss 
Percent  energy  loss 
Point  mass  contacts 
Point  mass  comes  to 
Point  that  force  is 


spring 

rest 

zero 
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loss  I Kl  t(qi  ~ q ^ 2 


100 


(q  - q * ) 2 ] - \ K 

7 Kj  [(qi  - q ") 2 - (q  - q')2] 

L o 


[ Cq i - q')2  - (q  - q')2] 

o 

7 Ki  [ Cq  l - q " ) 2 - q * q')2J 

L O 


or 

(3.1)  - K,  Cl  - ^T^-) 


Due  to  computer  limitations  we  must  approximate 
"F"  by  a continuous  function  and  where  there  is  jump 
discontinuities  we  use  a steep  linear  segment. 


where  e - small 
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Nominal  values  for  forward  buffer: 


FPLDCF) 

qo(F) 


* 2600  lbs 
= 5.2  in 


q 1 (F) 

% loss  £F j 
Kl(F) 

Ko(F) 


is  the  displacement  (dependent 
on  impact  velocity,  computed  in 
program) 

25% 

500  lbs/in 

computed  in  program  using 
equation  (3.1) 


C(F) 


* 0.01  in. 


fr3 


F 


(This  force  computed  by  a subroutined,  SPRG) 
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Branch  2: 


(3.1. a)  F - F ♦ (x  - x ) Kj 

O O 


Branch  1: 

Using  the  two-point  form  for  a straight  line 
with  the  end  points,  (x  , o)  and  (x  + e,  F + eKi) 

o o o 

(F  + eKi)  (x  - x ) 

(3 . 1 . b)  F = — 2 — 


Branch  5: 


, . „ , % loss 

let  G = 1 - 1Qo  ■ 


(F  G ♦ eGK ! ) (x  - x ) 
(3.1.C)  F * —2 


Branch  3: 


End  points  * (xm,  F_  + (xm  - x_)  Kj)  and 


m 


(xm  - c,  F G ♦ c*m  - c - X ) G Ki) 


(3.1.d)  F - F G * (x  - xm  * e)  [K,  (xm  - X ) *1^5'S  ♦ K,eG 


* r 1T&PJ  / • ♦ (x.  - xo  - .)  K,G 
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Branch  4: 


(3 . 1 . e)  F « F G ♦ (x  - x ) K G 

v O o 1 


4.  Force  from  rearward  buffer  (see  forward  buffer 
analysis . 

Nominal  values  for  rearward  buffer. 


FPLD(R) 

= 

2600  lbs 

qo(R) 

= 

0.2  in  (-  q is  the  displacemen 

qi(R) 

= 

dependent  on  impact  velocity 

* l0SS(R) 

X 

25% 

Kl(R) 

s 

500  lbs/in 

Ko(R) 

= 

computed  in  program  using 
equation  (3.1) 

C(R) 

s 

0.01  in 

fr4 

= 

F (This  force  is  computed  by 
subroutine,  SPRG) 

(6) 

fr5 

a 

- FSX 

fr6 

s 

- fs2 

(7) 

fr7 

s 

- fs3 

(8) 

fr8 

s 

- fs4 

(9) 

fr9 

s 

- FS  5 

(10) 

FRio 

3 

- fs6 

(ID 

FRu 

3 

- fs7 

(12) 

FR  1 2 

S 

- FCj 
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(13) 

(14) 

(15) 


FRi3 

fr14 

FRis 


- - FC2 

= * fc3 

= Piece-wise  constant  force 

acting  on  the  receiver  (accounts 
for  coulomb  type  friction) 

See  page  116 
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Forces  Acting  on  the  Slide: 

1.  Drive  Spring  Force,  FS i : 


q2 


Q 2 0 


Reference  position 


FS 


REF 


*SI 

Nominal  values: 
R 2 0 


FS 


REF 


VSI 


- Force  at  reference  position 

- Spring  rate 


= Calculated  (see  page  115  ) - 
Full  forward,  slide,  in 
equilibrium 

= 290  lbs 

**  13.2  lbs/in 


(Equal  and  opposite  forces  act  on  the  receiver,  see  FR5) 


> 
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- Distance  slide  travels  from  full  forward 
position  prior  to  gas  being  cut  off  to 
cylinder . 

6^  - Distance  traveled  (of  slide)  before  flow 

can  begin  to  be  drained  back  to  the  barrel. 


Pjj,  - Pressure  and  temperature  of 

the  barrel  gas. 

Pc.  Tc.  Pc  ' Pressure,  temperature  and  density 

of  the  cylinder  gas. 

P0,  T0*  PQ»  v0*  M0  ’ Pressure,  temperature,  density, 

velocity,  and  mach  number  at 
throat . 


- Throat  area  (orifice  area) 

- Orifice  coefficient 

- Flow  rate 

- Gas  constant 

- Ratio  of  specific  heats 
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(3.1.1) 


1 


Yb  ! 
' Rb  Tb' 


P M 
o o 


- 1 


M 2.** 
o i 


(See  "Compressible  Fluid  Flow,"  Shapino,  page 
for  equation  (3.1.1)<) 


(3.1.2) 


b 

F~ 


- 1 


- 1 


= (1 


Mq2) 


(3.1.3) 


M 2 
o 


Energy  equations  for  a deformable  control  volume 
(See  "Foundation  of  Fluid  Mechanics  and  Fluid  Statics 
Arthur  G.  Hansen,  page  130.). 


(3.2)  Q WSHAFT  + WSHEAR  + Jt  fVQ  v P (U  + y + gz)  dV 


+/s  p (h  + j-  + gz)  v ■ n dA 
c.v. 


c.v. 


p V, 


n dA 


Q 

w 

”SHAFT 

WSHEAR 

P 

U 


- Rate  of  heat  transfer  to  system 

- Shaft  work  being  done  by  system 

- Shear  work  being  done  by  system 

- Density 

- Internal  energy  per  unit  mass 
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V 

h 

gz 

n 

v 

r 

dA 

d 


- Velocity 

- Enthalpy 

- Gravitational  potential 

- Vector  normal  to  control  volume 

- Velocity  vector  relative  to 
control  volume 

- Velocity  vector  of  deformable 
volume 

- Element  of  area 

- Time  derivative 


. v. 


c .v. 


- Integral  over  control  volume 

- Integral  over  surface  of  control 
volume 


We  assume, 


(3.3) 


gz  " Q WSHAFT  = WSHEAR  ' 0 


V ^ 

Also,  y-  + gz  is  small  compared  to  U inside 
the  control  volume  (cylinder) . 


106 


(3.4) 


A 


Vc 

V 

c 

q 

c. 


'Pb 


- Initial  cylinder  volume 

- Cylinder  volume 

- Velocity  of  the  slide 

- Specific  heat  at  constant  volume 
of  barrel  gas 

- Specific  heat  at  constant  pressure 
of  barrel  gas 

- Internal  energy  per  unit  mass  of 
the  cylinder  gas 

- Cylinder  pressure 


The  energy  equation  is  given  by: 


d 

ar 


(ecpc  V 


W 


w 


o 


0 


(NOTE:  Ideal  gas  h - C^T) 
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or 


(3.5) 


dec 

3t~ 


P V ♦ 
c c 


e W 
c o 


VpJo 

b 


0 


Also , 

(3.6)  Pc  = (paV0  + /WQdt)  / Vc 


and 


(3.?)  Vc  = Vo  - ^ - <W  Ac 


Since  the  initial  number  of  moles  of  gas  in  the 
cylinder  is  so  small  compared  to  the  amount  of  gas 
that  flows  in,  we  can  assume  the  gas  constants 
(Cy,  Cp  and  R)  of  the  "initial  gas"  is  the  same  as 

the  barrel  gas  (gas  constants  in  cylinder  do  not 
change) . 


(3.8) 


(3.9) 


(3.10) 


Pc  • VJc  ■ <CP.  - cv  ) °c  Tc 

b b 


e = C T 
c vb  c 
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There  is  no  flow  if  the  slide  is  in  a given 
reg ion : 


If 

ab  1 - Cq2  - ^ref^  1 6a  » 


then 


0 


Position  of  slide  when  it  is  in 
equilibrium  position  (forward 
position  and  not  elevated)  (See 
page  115  ) 


(3.5)  along  with  equations  of 
motion  and  the  equations  that  describe  the  thermo- 
dynamic properties  of  the  gas  in  the  cylinder  determines 
the  system  dynamics. 

(3.13)  FS3  = (Pc  - Pa)  Ac 

at 


FS 


qREF 


REF 
*1  (SF) 


Energy  equation 


(3.14) 

(3.15) 

(3.16) 


T - T 
c a 


P , P 

c a 


t = 0 

atmospheric  temperature 
atmospheric  pressure 

Pa  1 (R  T ) 
c 


no 


Equal  and  opposite  force  acts  on  the  receiver 
(See  FR7). 

Nominal  values: 

Ac  - J CO . 588) 2 = 0.271  in2 

6 = 0.17  in 

a 

6^  = 3.625  in 

Ao  = (0 . 12) 2 = 0.033929  in2 

V - 11  J CO  - 588)  2 + 1.62  ~ (0.847)  2 + 0.375  (1 . 218) 2 
= 4.337  in3 

cPb  ■ °-4326  ■ 2-811  x 106  (H?)2  k 

Yb  - 1.31 

V • cpk  / *b  ■ 2-1459  x 106 

D D 

Pa  = 14.7  psi 

T = 530  °R 
a 
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Barrel  Pressure  and  Temperature  Determination: 


(3.17) 


We  use  an  experimental  pressure  curve.  Temperature 
during  "in-bore"  time  is  obtained  from  interior  ballis- 
tics calculations  (see  results  from  Darrel  Thompson 
and  Bill  Leech).  Once  the  projectile  leaves  the  barrel, 
and  assuming  no  heat  loss  and  using  the  fact  the  gas 
does  no  work  from  this  point,  we  see  the  adiabatic 
expansion  law  holds.  We  compute  the  temperature  of 
the  gas  from  the  experimental  pressure  curve  and  the 
following . 


of  gas  when  projectile  leaves 


gas  when  projectile  leaves 


Gas  begins  to  flow  after  the  projectile  reaches 
the  port. 


^b*1 


b_ 

bm 


bm 


bm 


Temperature 
barrel . 


bm 


Pressure  of 
barrel . 


tpoRT  = 0.00942  - time  to  reach  port 
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Barrel  Stagnation  Pressure  and  Temperature 


* 


Time  (Sec) 

Pressure 

Temperature 

0 

14.71 

530 

0.0001 

23300 

3741 

0.00062 

1600 

3411 

0.001042 

11500 

3323 

0.00118 

10510 

EXIT 

0. 00268 

5920 

0.00418 

3740 

0. 00S68 

2800 

ADIABATIC 

0.00718 

1765 

0.01168 

831 

0.01768 

520 

0.02068 

364 

0.050 

14.71 

Pbm  = 11500  Tbm  " 3325 

4.  Slide  rear  buffer  force  - FS4 


See  pages 
Nominal  Values 

95 

- 98  for  the  type  of  spring  force 

FPLD(SR) 

s 

0.0  lbs  (guessed) 

qo(SR) 

m 

9.817  in  (-q  is  displacement) 

* 

qi (SR) 

- 

(Dependent  on  impact  velocity) 

% loss 

s 

25% 

> 

Kl  (SR) 

s 

64000  lbs/in  (guessed) 

C (SR) 

s 

0.01  in 

fs4 

■ 

F 

(Force  computed  by  subroutine,  SPRG.) 
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Equal  and  opposite  force  on  the  receiver  (see  FR8) . 


5.  To  allow  for  time  for  the  round  to  be  fed  into 
position  in  front  of  the  chamber,  the  chamber  becomes 
seared  and  sear  is  released  when  the  receiver  is  some 
distance  rearward  in  recoil.  The  buffering  of  the  slide 
when  it  is  restrained  from  moving  forward  during  feeding 
is  accomplished  by  part  of  the  same  buffering  spring 
that  stops  the  slide.  (See  pages  96  - 97  for  type  of 
spring.) 

Nominal  Values: 


FPLD(SS) 

— 

0.0  lbs 

qo(SS) 

= 

9.817  in  (q  is  the  displacement 

2 

qi(SS) 

= 

(dependent  on  impact  velocity) 

Kl(SS) 

= 

48000  lbs/in  (guessed) 

% loss 

= 

25% 

e(SS) 

s 

0.01  in 

FS 

= 

-F 

5 


Equal  and  opposite  force  on  receiver,  FR^. 

FS  - F if, 

5 

qj  - q ( ] ) RELEASE 

q2  - q ( 2 ) LATCH  (during  cycle) 

q(j)RELEASE  - distance  receiver  is  rearward 
before  slide  is  released. 

q , -v  iatch  " distance  rearward  the  slide  must 
be  before  it  can  become  seared. 
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6.  Forward  Buffering  Spring  - FS 

6 

When  the  weapon  is  fully  locked,  this  spring  force 
restrains  the  slide  from  moving  further  forward  rela- 
tive to  the  receiver.  To  start  the  system  integration, 
this  spring  must  counter  balance  the  drive  spring 
force  along  with  the  weight  component. 

Nominal  Values  : 


FPLD(SF) 

— 

0.0  lbs 

qo(SF) 

= 

0.0  in  (q  is 

2 

the  displacement) 

qi(SF) 

= 

(dependent  on 

impact  velocity) 

Ki(SF) 

= 

64000  lbs/in 

e(SF) 

= 

0.01  in 

% loss 

= 

25% 

FS 

6 

= 

-F 

(Equal  and  opposite  force  acts 

on  receiver,  FR  ) 
1 0 

We  now  compute  q 
ZSL 


So  = (FSREF  - Ms  8 Sin  eo)/K1(SF) 


where 

e - Initial  angle  of  elevation 

. 0 

FSref  - See  page  102 

This  value  for  q insures  no  unbalanced  internal 
» forces,  initially.  20 

NOTE:  To  calculate  q we  assumed  the  preload  on 

2 0 

the  drive  spring  would  not  change.  Since  the  forward 
buffer  spring  is  much  stiffer  than  the  drive  spring, 
then  the  penetration  into  the  forward  buffer  will  not 
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be  significant  enough  to  change  the  preload. 


If  we  fire  with  slide  rearward  we  assume  the 
slide  sear  is  released  when  receiver  sear  is 
released. 

7.  To  allow  us  to  look  at  the  effect  of  forces  due 

to  stripping  of  rounds  or  coulomb  friction,  we  define 

a piecewise  constant  force  that  acts  between  points 

q and  q f . . , % and  opposes  the  motion  of  the  slide. 
2I.1J  2 U T r j 

If  q < 0 

2 

FS7  = FFRICT 


i = 1 , 2 , . . . , 2n 
j = 1,  2,  . . . , n 
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If  q >0 

2 

fs7  = ' ffrict 

i = 2n  + 1 , . . . , 2n  + 2m 

j = n + 1 , - . . , m 


Associated  with  each  force  will  be  an  equal  and 
opposite  force  acting  on  either  the  receiver  or  turret. 
A control  number  is  associated  with  each  piecewise 
constant  segment. 

If  1 then, 

FR  = -FS 

I 1 7 

If  0 then, 

FR  =0 

I I 

(Nominal  value  of  force  is  zero.) 

(Total  number  of  segments  allowed  is  3.) 
n m <_  3 

Forces  Acting  on  the  Chamber: 

1.  Rear  buffer  spring  force  - FC^ 

(See  pages  96  - 97  for  the  type  of  spring  force.) 
Nominal  Values: 

FPLD(CR)  = 1700  lbs 

q = 7.312  in  (-qc  is  the  displacement) 

i loss  = 251 
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i (CR) 
(CR) 


3400  lbs/in 
0.01  in 


FC 


= F 


(Force  computed  by  subroutine,  SPRG) 


FR  ) 
1 2 


(Equal  and  opposite  force  acts  on  the  receiver, 


2.  To  account  for  viscous  type  friction  we  define 


FC 


-C  q 
cHc 


Nominal  Values: 


FR  ) 

] 3 


(Equal  and  opposite  force  acts  on  receiver, 


3.  We  define  a piecewise  constant  force  (similar  to 
page  116)  acting  on  the  chamber  (FC  ). 

3 


(An  equal  and  opposite  force  acts  on  the  receiver, 
FR  ) 

1 *4 

Nominal  Value  is  zero 


n ♦ m < 3 


4.  The  force  from  the  propellant  effectively  acts 
on  the  receiver,  since  when  the  weapon  is  fired  the 
chamber  is  locked  to  the  receiver. 


118 


I 


c(t) 

tFIRE 

tDELAY 


Chamber  pressure  force 
Time  of  fire 
Ignition  delay 


FC  = -c(t) 


Nominal  Values : 

tFIRF: 

tDELAY 

t 

l 


i.See  control  equation) 
0 . 0 msec 
10.0  msec 
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Effective  time 


t DELAY  ( 


0 

0.00075 

0.0015 

0.003 

0.0045 

0.006 

0.00675 

0.0075 

0.00767 

0.00787 

0.00802 

0.00839 

0.00895 

0.00919 

0.00942 

0.00984 

0.0105 

0.0120 

0.0135 

0.015 

0.0165 

0.021 

0.027 

0.030 

0.0305 


Force 
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Torsional  Forces  Acting  on  the  Chamber 


Other  than  constraint  torsional  forces  on  the 
chamber,  there  will  also  be  frictional  forces.  We 
define  a piecewise  constant  torsional  force  (similar 
to  the  translational  force  on  page  116)  acting  on 
the  chamber. 


piecewise  constant 


(Nominal  is  zero;  three  piecewise  constant  segments) 

Since  we  have  no  degree  of  freedom  for  roll,  the 
equal  and  opposite  torsional  force  applied  to  the 
chamber  by  the  receiver  does  not  need  to  be  considered. 

The  forces  acting  on  the  receiver  are  given  by: 


l FR- 


R 


Forces  acting  on  the  slide  are  given  by: 


Fy  = l FS. 
AS  j J 


Forces  acting  on  the  chamber  are  given  by 


FY  - E FC, 
AC  k k 


M 


= M. 


Cl 


The  only  thing  left  to  determine  at  this  point 
is  the  control  equations  that  determine  when  to  fire 
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Preliminary  Matrix  Relationships 


Part  A: 

If  A is  a column  matrix  and  a function  of  q (a 
column  matrix) 


Aj  (q) 

A (q) 

2 


An  (q) 


then  from  the  chain  rule  for  differentiation  of  a 
column  matrix 


dA  3A  da 
3t  = 7=  3t 


or  we  can  write  this  in  the  form; 
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t 


9Ai (q ) 
9qn 


( 


(1.2) 


dA 

Ji 


9 A ] (q)  9 A i (q ) 

3Qi  3q2 


9An (g)  9An (q ) 

9q i 3q2 


9An  (q ) 
9qn 


If  A is  a function  of  t,  q and  q then  we  can  see 
that , 


M 31  dA(t  .q,q)  . 9A  d£g.  + 9A  dg.  + 9A 
dt  dt2  9q  dt  9 1 

Part  B: 

Lagrange's  Equation: 


)> 


i = 1 , . . . n 

T - Kinetic  Energy 

V - Potential  Energy 

P - Dissipative  Potential 

Q - Generalized  Forces  (i  = 1,  2,  . . . n) 

qi 

q ^ - Generalized  Coordinates  (i  = 1,  2,  . . . n) 
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Using  notation  from  Part  A and  letting 


Lagrange's  equation  (equation  2)  becomes  in  matrix 
form: 


(2.1)  at  (1t)T  ' (i=)T  + 


(iI)T  + 

3q 


3q 


Qq 


To  integrate  the  equations  of  motion,  one  needs  to 
separate  the  second  order  derivatives  from  the  first 
order  terms.  If  we  apply  equation  (1.3)  to  the  term 


in 

is 


d f3T, 

at  (“) 


3q 

equation  (2.1)  (Note:  Assuming  that  kinetic  energy 

a function  of  t,  q and  q)  Larange's  equation  becomes 
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4 (4)T  q » 4 c4)T  k * 4 <4)T  - (4)T 


3 q 3q 


3q  3q 


T T 

- C^)  - c^r)  - Qq 


With  this  form  of  the  equation,  if  we  are  able 
to  determine  T,  V,  and  P as  functions  of  the  generalized 
coordinates,  we  can  use  an  "algebraic  computer  language" 
to  generate  the  equations  of  motion.  The  terms 

4 (1I)T  , 4 (il)T  , 4 (4/  , (4)T  , (4)T  . and 
3q  3q  3q  3q  3q  3q  3q 


are  easily  generated,  punched  on  cards,  or  saved  in 
some  other  form.  These  expressions  can  be  used  in  an 
integration  routine. 

The  term, 

— (— )T 
3 q 3q 

can  be  considered  a generalized  mass  matrix.  In  order 
to  compute  the  second  order  derivatives,  the  inverse 
of  this  term  must  exist.  We  denote  the  mass  matrix  by: 


M = — (— ) 

• ' • 

aq  3q 
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We  now  specialize  these  equations  to  a system  of 
rigid  bodies.  We  show  that  if  the  center  of  gravity, 
C.G.,  is  determined  relative  to  the  inertial  reference 
frame  in  terms  of  the  generalized  coordinates  and  the 
angular  velocity  of  each  body  relative  to  a fixed 
body  axis  and  in  terms  of  the  generalized  coordinates 
is  determined,  then  we  show  how  to  use  the  computer 
to  operate  on  these  expressions  to  generate  the 
equations  of  motion. 
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Inertial  coordinate  system 


X Y Z 

o o 

< 

X-Y.Z-  Body  axis  of  the  "i-th"  mass 

(Origin  at  the  C.G.  of  the  "i-th" 
. mass) . 


r^  Vector,  locating  C.G.  of  the 

"i-th"  mass  relative  to  the 
inertial  reference  frame. 

Angular  velocity  of  the  "i-th" 
mass  relative  to  the  body  axis. 


The  kinetic  energy  of  the  "i-th"  mass  can  be 
written  as: 


T T 

Ti  * 7 Mi  F1  ?i  * 7 wi 


I.  w. 
1 1 


where 

Mass  of  the  "i-th"  body 
Inertial  tensor  of  "i-th"  mass 


xixi  xiyi 

xizi 

i 

I 

Yixi  yiyi 

yi  zi 

I 

I 

zixi  ziyi 

zizi 

(Note:  li  = I.T) 
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Identify  and  as  being  column  matrices: 


also  and  are  functions  of  q,  q and  t: 


r = r (q,  t) 

^ = W (q,  q,  t) 


The  total  kinetic  energy  is  given  by: 


. T • T 

.3)  T - I Tt  - I (}Mj  Fj  rt  • 1 I4  W.) 


Also,  r^  is  given  by: 


3r . 
l 

3t 
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We  can  now  rewrite  equation  (2.2)  as: 


3 3T-  T 1 9T . T • 3T- 

(2.4)  (-4)  q + * -4  (-4)  q + £ — ~ (—4) 

i i , — i 3t 

3q  3q  3q  3 q 3q 


3T- 


? (—4)  + (4)  + c-4)  = Qq 

3q  3q  3q 


The  generalized  mass  matrix  is  given  by: 


Z 

i 


M. 


l 


1 3q  3q 


Lower  order  terms  involving  are  given  by: 

T T T 

, 3T,  • - 3T.  3T. 

r C.  - E [—2-  (-4-)  q + aT  (— t)  ■ (— ) j 

1 1 3q  3q  3q  3q 


The  mass  matrix  of  the  system  is  just  the  sum  of 
mass  matrices  from  each  body.  Also,  the  first  order 
terms  are  given  by  the  sum  of  the  lower  order  terms. 

We  can  derive  equations  of  motion  b^  looking  at 
a particular  mass  and  determine  r.  and  w.;  use  the 
computer  to  determine  M.  and  . 1 1 

(NT  - mass  matrix;  - column  matrix  of  lower  order 

terms . ) 
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These  expressions  are  printed  out  on  cards  or 
some  other  convenient  form  and  used  in  an  integration 
routine.  During  integration  of  equations  of  motion 
by  numerical  methods  we  compute 


(2.5)  q = M'1  (-  iil  - — - IC-  + Q— ) = f (q,  q,  t) 

— — i ^ 

3q  3q 

If  we  use  a routine  that  solves  a system  of 
the  form; 

y = fi  (x,  y) 

(2.6)  y (x  ) = y (Initial  conditions) 

O O 

then  we  need  to  transform  (2.5)  into  that  form. 
Let 

(2.7)  q i - q 

S = ^ 

(2.8)  y = 
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(2.9) 


( 


r 

(2.10) 


(2.11) 

1 


1 


y 


(y,  t) 


(Initial  conditions) 


So  far  nothing  has  been  said  about  generalized 
forces.  If 


is  the  point  of  application  of  a force 


then , 


9r 


Qq 


z 

p 


f - 

_£  F 


aq  p 


(Note:  From  literature  Q 


) 
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We  can  use  the  computer  to  derive  the  equations 
that  determine  the  generalized  forces.  Read  in 

?£■  and  F , and  compute  Q-  . 

P P 


We  now  obtain  the  expression  for  the  generalized 
force  when  a moment  is  applied  in  the  direction  of  a 
given  axis. 


We  will  show  that, 


T 


30 

3qi 


Replace,  torque  t by  the  couple  (see  diagram  below). 
t = 2 | | F | | II  rs  ||  Sin  0 

t - 2 ||  F ||  ||  Fs  || 
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Since  F is 


and , 


(2.13) 


(2.14) 


(2.15) 


(2.16) 


II 


II  " 


II 


r 


s 


II 


then 


F I I 


From  equations  (2.12)  and  (2.14)  and  the  fact 
x = 2 | | F ||  | | rs  || 


^6  is  a function  of  the  generalized  coordinates 
and  the  computer  can  be  used  to  generate  (2.15). 

In  vector  notation  (2.15)  is  written  as: 


T 


36 

3q 


134 


